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The large theory is constructed for the metalhc state of high-T^. cuprates. The Emery three-band model, 
extended with 0,-Ov hopping tpp, and with Ud — » is mapped on slave fermions. The Dyson time-dependent 
diagrammatic theory in terms of the Cu-0 hopping tpj, which starts from the nondegenerate unperturbed ground 
state, is translationally and asymptotically locally gauge invariant. However, it replaces the anticommutation of 
the fermions on the Cu and O sites by commutation, so it is antisymmetrized a posteriori. The small parameter 
of the theory is nj. The lowest order of the theory generates the single particle propagators of the hybridized 
pdp- and rfpd-fermions with the exact covalent three band structure. The leading many-body effect is band nar- 
rowing, obtained without mean-field slave particle approximations. It is accompanied by Landau-like damping 
of the single particle propagation, due to incoherent CU-O2 intracell "mixed valence" fluctuations. The corre- 
sponding continuum falls below the Fermi level, which lies in the lowest hole band. The conventional Luttinger 
sum rule for coherent bands is thus broken. Due to local gauge invariance, these results are insensitive to the 
omission of the Cu-0 anticommutation rules. The latter affect the effective local kinematic repulsion Ud/i be- 
tween the hybridized pdp hole propagators. Its a posteriori antisymmetrization removes the repulsion between 
the pdp particles in the triplet configuration, keeping their singlet repulsion intact. Such an Ud,, is the metallic 
counterpart of the superexchange Jpd- Udp is weak to moderate, and close to optimal doping favors the coher- 
ent, antinodal, low energy SDW correlations, remotely related to the formation of the marginal Fermi liquid. 
Resonant valence bonds which involve Jpd appear as incoherent perturbative corrections to these correlations. 
In function of doping and/or frequency, the incoherent "mixed valence" fluctuations compete with magnetic co- 
herence. This dichotomy, between properties local in the direct and reciprocal space, is the fundamental feature 
of the Emery model in the large Ud metallic limit. 



I. INTRODUCTION 

The long standing question in high-Tc cuprates concerns 
the nature of interactions which are responsible for supercon- 
ductivity and other unusual properties of these materials. The 
early high-energy spectroscopic measurements indicated that 
the Hubbard interaction Ud on the Cu site might be quite 
large.' This opened the question of whether or not Ud, to- 
gether with the appropriate set of the single particle param- 
eters, can account for correlations observed in the high-Tf. 
cuprates. A related question^ ^ is what are effective interac- 
tions which include the strong electronic correlations in high- 
Tc cuprates. The fundamental aspects of these questions are 
discussed here from the theoretical point of view. 

The observed phase diagram of cuprates is well known.Gl^It 
is characterized by the principal crossover (ignoring the small 
interplanar couplings) between the insulating long or short 
range AF-Mott-like phase at small hole doping < x < x^s 
and the metallic phase for x > Xcs- The latter shows coher- 
ent and incoherent features in single particle propagation and 
in Raman, optic and magnetic correlations^ The coher- 
ent single particle propagation emerges at x^s from the broad 
background, sharply on arcs around the nodal points in the re- 
ciprocal space, as best observed by ARPES." The metallic 
properties of cuprates stem from such nodal behavior Con- 
comitantly, the leading harmonic of the coherent magnetic 
correlations, incommensurate for x > 0, rotates at x^s from 
the diagonal to the collinear position with respect to the main 
axes.^This is reflected in the ARPES spectra for x > Xcs in 
the appearance of the pseudogap in the vicinity of the anti- 



nodal van Hove (vH) points, which are on the other hand 
inessential for metallization. Typical experimental values of 
Xes found from all these (and many other) measurements fall 
below five percent. These measurements provide the average 
hole occupations of the Cu and Ox,y sites rid and rip. I n the 
X > Xcs metallic phase, averag e site occupations founcP^EIl 
by NQR in cuprates corresponciSEll to rid appreciably less 
than unity, and increases on doping with holes. An ac- 
curate measurement of rid in the x < x^s regime is difficult 
in the AF phase, but it is usually believed'^ that in the range 
< X < Xcs rid decreases slightly on doping. Concomitantly, 
the chemical potential observed by XPS decreases smoothly 
with X through the crossover at x ^ XcjI^^^The early soft- 
X-ray measurements on lanthanates^"^ indicate that the smooth 
behavior of the chemical potential is accompanied by the pro- 
gressive transfer of the spectral weight to lower electron ener- 
gies. 

In view of these and other experimental results, particu- 
lar care is devoted here to the identification of the theoretical 
regime which is appropriate for cuprates. Two main doping 
regimes are distinguished with respect to Xcs- The present pa- 
per is focused on the x > x^., metallic regime. The simultane- 
ous appearance of the single particle band structures together 
with magnetic correlations which tend to rotate in "collinear" 
Bragg points is a clear signature of the charge and spin co- 
herence in the low energy eigenstates of the x > x^j metallic 
regime. This regime is usually further divided in underdoped, 
optimally doped and overdo ped ran ges. The superconductiv- 
ity, which occurs exclusively ! ' 3 | 25 | 26 | ^^j. ^ ^ ^^^^ itself a co- 
herent state. Moreover, these coherent features, "localized" 
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in the reciprocal space, are not only balanced among them- 
selves but also compete with dynamic charge and spin disor- 
ders within the CuOa unit cell, localized in the direct space. 
The main aim of the present work is therefore to improve the 
theoretical understanding of such competitions. 

The presentation of our work, reported preliminarily in 
Ref '2T is divided in two papers. After comparing the weak 
and strong coupling limits, the large perturbation theory 
is developed for x > Xcs in the present paper The theory is 
formulated at the outset as an expansion procedure in terms 
of the small average occupation n^' in the Hartree-Fock (HF) 
3-band state. The companion paper compares these and ad- 
ditional theoretical results to the experimental data. Even- 
tually it will be argued there that the results derived in the 
present paper agree qualitatively and often semi quantitatively 
with the measurements in the well developed metallic phase 
of cuprates and especially in the range around and beyond op- 
timal doping. 

II. WEAK VERSUS STRONG COUPLING 

A. Theoretical model 

The theories of high-T^, cuprates often start from the tight 
binding model, with the vacuum consisting of Cu^(d'") and 
(P^) states. The Cu^^(d') and 0"(p^) site energies are de- 
noted then by and Sp with Apj - Sp - > in the hole 
language. The Cu-'^(d^) state is reached by spending the en- 
ergy 2sd + Ud where Ud describes the bare interaction of two 
holes on the Cu site, which is reduced to some extenPby in- 
tra atomic correlations. Referring to LDA results, the ©(p"*) 
configuration is usually associated^** with the energy 2ep i.e. 
Up is considered as relatively small. This was completed 
with the Cu-O hybridization tpd and often called the Emery 
model.^"* The original model was later extended^" by the di- 
rect O-O hopping tpp which describes the hole propagation 
rotated by njA with respect to the CuOa axes. The model is 
completed by fixing the total number of holes \ + x per CuOa 
unit cell, where x is the number of doped charge carriers, as- 
suming that X < 1 . The average single particle occupations of 
the Cu and Ox,, sites and Hp are then linked by the sum rule 
Hd + 2np - \ + X, i.e. - lup or «,/ itself is the "primary or- 
der parameter" of cuprates understood as CU-O2 charge trans- 
fer (CT) salts. Together with the associated "mixed valence" 
fluctuations, «^ will play a prominent role here. 

B. Properties of weak coupling results 

It will appear below that under appropriate conditions the 
coherent features of cuprates can be explained within an ef- 
fective weak coupling theory. This was first attempted"*^ un- 
der the assumption tpd < Ud < Apd, using an effective in- 
traband Uj reduced by metallic kinematic^"* correlations. It 
is therefore interesting to review briefly the properties of the 
small Ud theory. Let us first assume^ that tpp = 0. At T=0 
the X = Fermi surface passes then through the logarithmic 



vH singularities, irrespective of the ratio tpd/Apj. The equal 
sharing of charge = 1/2 at x - b etween Cu and two 
O's is obtained in the covalent 

limipm 

tpd » Apd, whereas 

Hd - \ corresponds to the opposite ionic limit Apj » tpd- The 
X - Q Fermi surface is perfectly nested and this may lead to 
marginal Fermi Ik^uid.'^^ Finite Ud, and especially its Umk- 
lapp component,^"* enhances quite strongly the antiferromag- 
netic (AF) correlations with commensurate q^p - Gjl SDW. 
The ^ = O JOy CT fluctuations within the CUO2 unit cell are 
also enhanced. These fluctuations are coupled quadraticalljl^H 
to the very slow tilting modes in lanthanum cuprates which 
may resulPEili ^Yie LTO/LTT instability. With finite dop- 
ing X the SDW instability moves^** ^'^ to incommensurate val- 
ues of qs Dw and becomes weaker. In particular the effect of 
the G-Umklapp interaction Ud in the build up of AF correla- 
tions is diminished in this way. Usually, the Umklapp in ques- 
tion is removed by hand from the theory when qs dw becomes 
incommensurate, i.e. the SDW commensurability pinning is 
ignored. This results in a (too) smooth sliding of qs dw on in- 
creasing X. Although the SDW behavior emerges more or less 
correctly, the main problem of this description is that it putJ^H 
the commensurate LTO/LTT instability at x - Q rath er than at 
sizeable hole doping x ^ x,,// ^ 1/8, observed^^^ED in LBCO 
when the Fermi level reaches'*^ the vH singularities. 

This could be remedied by including tpp?'^ Although 
smaller than tpd on chemical grounds, it is quite rele vant in the 
weak coupling theory. Fir st, for ?pp < (appropriatd^^EIllIlfor 
high Tc cuprates), tpp sets ■'^'^IIllEI the Fermi level of the x = 
half filled lowest band below the vH singularity, which means 
that the latter is reached upon a finite hole doping Xyfj > 0. 
The ARPES spectra of the ho le-do ped cuprates in the x > x^s 
metallic state can than be fil^EIl by the Emery three-band 
structure. These fits indicate"*^ that the bare band parameters 
obey the relation A^^^ > 2?^^ > Apd\tpp\ with \tpp\ large enough 

to account-" for the n/A rotation of the Fermi surface (Fermi 
arcs) with respect to its tpp - form. Concomitantly, the Cu 
occupation «j at x = is reduced below 1 /2 in the covalent 
limit tpd > Apd- Furthermore, tpp breaks'*"*'*^ the perfect nest- 
ing properties of the X = XyH Fermi surface, i.e., it plays the 
role^^ES) of the imperfect nesting parameter. The elementary 
SDW particle-hole bubble develops then the peaks at incom- 
mensurate qs DW for X - XyH and small (ypnEllEll When the 
small interaction Ud is introduced, the resulting G-Umklapp 
scattering of two particles is in discord with this value of the 
wave vector and the resulting value of qsow = (^0 + G)/2 
is generally incommensurate and weakly affected by a small 
Ud- qo characterizes also the lattice deformations (stripes^ 
which accompanies the SDW.^^ In parallel, the elementary 
O JOy CT particle-hole bubble is^^ (logarithmically) singular 
for X - XyH at w = 0, ^ = for a ny val ue of tpp. This cor- 
responds to the Jahn-Teller splitting'^^'^ of the vH singulari- 
ties, and presumably favors the commensurate LTT instability 
in some lanthanates. All mentioned features agree well with 
observations^ . In particular, the results of the single band 
t, t', t", Ud weak coupling model obtained in the improved 
ladder approximations'** are widely used in fitting the ob- 
served magnetic features of the metallic phase in cuprates. 
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However, the problem is that strong magnetic correlations 
occur in the weak coupling theory only for x ^ x^h- As men- 
tioned above, this particular property is at variance with ob- 
servations that appreciable magnetic correlations (coupled to 
the lattice, i.e. stripes**) persist over a wide range of doping, 
from jc > up to optimal doping and beyond. While the weak 
coupling theory which uses finite tpp thus explains the mag- 
netism of the well developed metallic phase reasonably well, 
it fails to describe the emergence of the Mott-AF phase on 
decreasing the doping towards x < Xcs- 

Another fundamental feature not explained by the weak 
coupling theories is the appearance for x > Xcs of the broad 
background in the single particle propagation and related cor- 
relations. Not to be confused with the marginal Fermi liquicP^, 
this background is associated here to the intrinsic intracell, 
incoherent "mixed valence" d"'-i-p^«->d^-i-p'' charge transfer 
(CT) fluctuations which involve, bes ide the two O-sites, the 
Cu-site with large interaction f/j.UJ^This generalizes the re- 
sults for the Anderson lattice where the mixed valence fluctu- 
ations involve the site rather than the unit cell. Already in the 
earliest theoretical works on cuprates with Ud » Apd, it was 
pointed out"*** that the static d'"«->d'' disorder is an essential 
feature of the tpd = limit when sj falls within the disper- 
sive band, Ap^ < 4\tpp\ in the present language. Actually, for 
large Ud the finite tpd is bound to render this local disorder 
dynamic^^. Indeed, the broadened single particle propagators, 
named here pseudoparticles (in distinction to quasiparticles) 
by analogy with the Hubbard model,i^, appeared in the early ^' 
large Ud slave boson NCA calculation, in high temperature 
expansions'^ and also in the non-magnetic'^' DMFT calcula- 
tions with Ud > Apd- However, the direct relationship between 
the broadening of the single particle spectra obtained in these 
calculations and the dynamic d'"«-»d^ disorder of the large Ud 
limit was recognized only recently and is further developed 
here by emphasizing its connection to the Landau damping by 
local (intracell) "mixed valence" fluctuations among the Cu 
and O sites. 



C. Strong coupling limit 

1. Mean-field slave boson theory 

The described mixing of weak- and strong- coupling fea- 
tures motivates us to investigate carefully the Ud » Apd 
limit of the f^^ -extended Emery model. In our endeavors 
we are led to some extent by the translationally invariant (i.e. 

"metalh c") U d - mean field slave boson (MFSB) theory,'"* 
extende(p55l 

with harmonic fluctuations of the boson field 
(MFSB-nfl). It is therefore clarifying to survey briefly these 

results here. 

The MFSB with tpp - was first usecPEH to describe 
at jc = the change of an insulator into a correlated metal 
through the Brinkman-Rice (BR) "phase transition" between 
the insulating n'f^'^* = 1 and metallic < n'f^'^* < 1 
states.lS^. For \tpp\ smafl with respect to t^pd/^pd the BR tran- 
sition is shifted linearly in tpp to higher values of (A^^/f^j) 
and for large \tpp\ it moves to 4-\tpp\ » Apd^^ This is illustrated 




Figure 1 : (Color online) Position of the MFSB x = transition be- 
tween metal (yellow) and BR insulator (green) in the tpp/tpj, Apj/tpj 
parameter space. The insert shows schematically a characteristic be- 
havior of 1 - n^''^^ at X = in function of Ap^/tpj for a given small 
tpp/tpd (yellow and green). The smooth behavior expected here in 
absence of magnetic correlations is also shown (black). 

in Fig. [T] together with the corresponding behavior of the CT 
order parameter n^*^^'^*' atx = 0. n'"^^^^^ is intimately related 
to the renormalized hopping f^^^(l - n^*^^'^^'). For finite hole 
or electron doping x, the MFSB smoothes out the x = BR 
transition. Above the transition a small doping of the x = 
state produces dn'^^^^ Idx < and below it dn^''^^/dx > 0. 
At finite x the change of sign of dn'^^^^ jdx occurs below the 
BR point.'* Noteworthy, n^'^^^/dx - corresponds to the 
situation when the additional hole goes at average to the oxy- 
gen sites without leaving the renormalized lowest hybridized 
band. 

Although the BR transition is an artifact of the mean-field 
approximation this notion will prove useful in distinguishing 
the three parametric regimes of interest. Actually, the regimes 
"well below", "well above" and, transitory, "around" the BR 
transition will remain distinct in the present approach. 

Let us turn next to the MFSB-i-fl. The local gauge in- 
variance, although satis fied o n average,'^ " is irremediably 
broken in the MFSB-nfl.EHl^ The deep flaw of the MFSB, 
MFSB+fl (and to lesser extent of the early slave boson 
NCaP^SI) is that they start from the unperturbed HF state of 
extended spin carrying slave fermions, rather than from local- 
ized ones. This choice breaks the exact local gauge invari- 
ance at the outset. Moreover, through /g = 1 - n^'^^^ the 
MFSB and MFSB-nfl allow for finite average boson displace- 
ment /() - {f + f^) which is manifestly forbidden by the local 
gauge invariance. In addition, the MFSB+fl introduces the 
harmonic fluctuations of the boson field asymmetrically with 
these of the fermion field. The "mixed valence" fluctuations 
are thus poorly approximated in the whole range of parame- 
ters. 

Well above the BR transition the MFSB-i-fl produces the 
narrow HF "resonant band" with a reduced projection"*^ " on 
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copper site. This band contains 2x states which accommodate 
X additional holes essentially on oxygens i.e. close to energies 
Bp. Such band is thus approximately half-filled, i.e. pinned to 
the Fermi level. The rest of the carriers are in the dispersion- 
less, incoherent states close to e,/. 

In contrast to that, well below the BR transition, the 
MFSB-Hfl renormalizations become weak. The renormalized 
band parameters t^f^^ and t^'^J^^ are weakly dependent on 
X, while Bp and tpp are 

unaffected.Esmil] 

The Fermi energy 

lies, after phenomenological correction,"*^ above the remnant 
incoherent contribution below and the band accommodates 
therefore somewhat less than \ + x carriers. This looks very 
much like the usual band theory with weak magnetic corre- 
lations neglected, in agreement with expectations that despite 
large Ud the renormalizations are weak below the BR tran- 
sition. However, the average boson displacement /o is large 
(close to unity) below the BR transition. The corresponding 
renormalization effects of the HF state by the MFSB-nfl, al- 
though small, are therefore suspicious. 

The problems of the MFSB-i-fl become more visible on ap- 
proaching the BR transition. The harmonic fluctuations of the 
boson field do not aff'ect the condition for average local gauge 
invariance, because the associated number of shifted spinless 
bosons vanishes. The MFSB+fl preserves thus the BR transi- 
tion instead of smoothening it out, as could be expected from 
the "mixed valence" CT fluctuations conjugated to the CT or- 
der parameter 

This is particularly unsatisfactory because the behavior in 
the vicinity of the BR transition is believed to be relevant 
for cuprates.'I^'SlSl] Such assignment of cuprates relies mostly 
on the NQR analysi^^ leading to considerably less than 
unity with dudjdx > for x small. This behavior of n^ix) is 
consistent with recent measurementdS^Ei of d/u/dx, which is 
small throughout the crossover at small Xcs- The correspond- 
ing regime of bare parameters A^^ > 2f^^ > Apd\tpp\ > Atjp 
already mentioned at the beginning of this Section falls into 
the vicinity of the BR point and describes well the ARPES 
measurements in high-T^ cuprates. ^^1^2111 

Our endeavor can thus be understood as an attempt to red- 
erive results similar to these of the MFSB+fl by respecting 
better the local gauge invariance on improving the description 
of the "mixed valence" fluctuations and on including the miss- 
ing magnetic correlations. In the present paper this will be at- 
tempted on approaching the BR transition from below. It will 
turn out that "mixed valence" fluctuations, entirely neglected 
in the MFSB and inconsistently treated in MFSB-i-fl, play an 
important role in the perturbation theory,^^ especially close to 
the BR point. In particular, the BR "transition" will be re- 
moved by incoherent local "mixed valence" fluctuations. The 
BR collapse of the coherent band-width is thus absent. The 
detailed comparison of the (nonmagnetic) perturbation theory 
and the MFSB+fl is given in Sec. [Ill] 

Let us finally point out that our approach from below the 
BR transition bears resemblance, before magnetic correlat ions 
are included, not only with MFSB or MFSB+PEna and 
NCA,^' as further discussed below, but also with large Ud 
DMFI^^ ^' and LDA+lP^ calculations for cuprates. 



2. Mott and metallic limits with magnetic correlations 

The magnetic correlations can also be conveniently dis- 
cussed in three regimes introduced above in Fig.[T] Well above 
the BR transition it is usual to start for small dopings from the 
unperturbed n™^ = 1 Neel ground state which corresponds to 
the average Cu occupation close to unity. That state is widely 
used to approach the propagation of the first additional hold^ 
or electron''' in the x=0 AF phase of cuprates. The d oped 
hole (i.e., the chemical potential for x = 0+) is placed^lsol 
on th e upper o xygen level close to Bp. Orig inally, the t — J 
modePEUmi with superexchange^iSillil y^^ = ^t^dl^ld 
was so obtained for Ud > t^pd tpd > and tpp = 0. Such 
Jpd is small well above the BR-point and in this sense it in- 
troduces small energy scale for AF correlations at small dop- 
ings. At large dopings the magnetic effects are progressively 
removed. 

In contrast, well below the BR transition in Fig. [T] f/j > 
tpd » Apd > and sizeable tpp hybridizations lead to a re- 
ducti on of rid to small values for x = 0. It is then preferable 
to puP3i3the first additional hole (i.e. //) on the lower (rather 
than upper) covalent level. Its average charge in the unit cell 
goes then mostly to the O sites. The single particle hybridiza- 
tions among the unit cells are further included on the equal 
footing. It is ensured in parallel that holes on O sites avoid 
each other on the intermittently empty Cu site in the tempo- 
ral dimension through their effective repulsion, as was em- 
phasized in the preliminary report of the present work.^^ This 
means that the t - J model^'' is not a suitable starting point 
for tpd s> Apd and sizeable tpp. Rather, the metallic phase, 
as described below through dynamic spin and charge trans- 
fer disorder associated with "mixed valence" fluctuations and 
magnetic correlations due to the effective repulsion^^ between 
holes on oxygens, prevails then even at x = 0. 



3. Intermediary regime 

The regime A^^ > 2f2^ > Apd\tpp\ > 4tlp, close to the BR 
transition in Fig.[T] requires additional care. The unperturbed 
Neel state is always appropriate around x > 0. At x = the 
unity of charge is however nearly equally shared, nd ~ 1 /2, 
between copper and two oxygens. The intracell covalent 
hybridizations,^^ bare tpd in particular, are thus included from 
the outset.'''^ The additional hole goes to the energies approxi- 
mately half-way between Bd and Bp » Bp - 4-\tpp\. The effects 
of superexchange between holes on coppers and effective re- 
pulsion between holes on oxygens compete further in func- 
tion of doping (and frequency). The approaches based on the 
magnetically ordered = 1 unperturbed ground states of 
the Ud - °° Emery model are then relevant at low frequen- 
cies (temperatures) and sufficiently small dopings. However, 
the crossover to the metallic state accompanied by the small 
variations of nd and p is expected^^ to occur for small x ^ Xcs- 

We therefore take that the approach based on the Neel un- 
perturbed ground state is restricted to the small range x < x^s- 
In the case when ju for x = 0+ is closer to Bp the Zhang-Rice 
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rij X If-/ approachP^mHU is a reasonable starting point 
at small x. However, the effective repulsion of the holes on 
oxygens via the intermittently empty copper site has to be in- 
cluded beside Jpj. In the opposite case the effective repulsion 
of the holes on oxygens might be a dominant interaction, as 
further discussed in the last Sec.lVIl 

On the other hand in the doping range x > x^s the distinc- 
tion between the two mentioned cases becomes less obvious. 
We will describe this range in terms of a renormalized three- 
band theory with "mixed valence" fluctuations and magnetic 
correlations, approaching Xcs from above. It will turn out [see 
in particular Sec. |VI[ that the coherent magnetic fluctuations 
due to the effective repulsion between holes on oxygens are 
enhanced with respect to the effects of Jpd by nesting of the 
vH singularities. The latter occurs in the range around optimal 
doping. Such behavior is better mapped o n the widely used f, 
f', t" model with weak effective repulsiorP^I^, than on the f, 
f', t" , J model. The former mapping results'^ in "coUinear" 
incommensurate magnetic correlations. Such metallic regime 
includes perturbatively the fluctuating valence bondings due 
to superexchange Jpj. On decreasing x towards x^s these cor- 
relations tend to the AP ordering of copper spins. The picture 
completed with magnetic correlations is discussed in Sec.[V] 



ni. SLAVE FERMION PERTURBATION THEORY 
A. General 

1. Hamiltonian and its Hilbert space 

It is well k nown th at the Uj - oo theory can be mapped on 
slave particle J5Sl23|74'^jj(j mapping is only briefly reviewed 

here. The d'" state on Cu at the position ^ is denoted by /!|0) 

R 

and the d'^ state with spin cr by ^|0), where |0) is the aux- 

a-R 

iliary vacuum on Cu. In the so spanned three-state space (d^ 
state with energy led + Uj is omitted at |x| < 1), the number 
operators of the slave particles satisfy Qg = n^g+Yjo- - 1- 

The physical fermion c' ^ projected on the d^, d'" subspace 

o-R 

is written as c'^ — » ^'^/r- The corresponding number op- 
erators satisfy n°"^ = n°"^, usually called the "Luttinger sum 

dR bR •' 

rule" (LSR). fe' ^ and can be taken respectively as fermions 
and bosons ("slave boson representation", SBR) or as bosons 
and spinless fermions ("slave fermion representation", SFR) 
in order to satisfy the anticommutation rules on and among 
Cu sites projected on the d^, d"* subspace. The states on oxy- 
gens are associated with physical p-fermions so that the SBR 
satisfies the Cu-O anticommutation rules while SFR requires 
additional care. Any physical Hamiltonian in the Ud - oo 
limit can be written in terms of p-fermions and these slave 
particles. It is locally gauge invariant and commutes with the 
operator Qj^, i.e. = 1 is a physical constant of motion. 

More explicitly, the slave particle Hamiltonian can be writ- 
ten as 



H, = Ho, - AN + Hi, 

HoA = Hodisd) + Bp, tpp) . 

(1) 

R 

When Ud - oo is the only interaction^^ Hj is proportional to 
tpd which originally describes the single-particle hopping. A 
enters the site energy of the /-and /j-particles appearing, as 
usual, by adding A{Qj^ - 1) into the unperturbed Hamiltonian 
Hq. Since such H(u is additive in three b-, f- and p- particles 
the associated full Hilbert space is the direct product of three 
single particle Hilbert spaces and the unperturbed ground state 
is 

|Go) = |G*)®|G^)®|Gj;). (2) 

According to Eq. (jTji Hodisd) + '^YiiiQi ~ 1) is diagonal 
in the site representation, i.e. the Cu-sites are either occupied 
by the /- or by the /j-particle (either d'" or d') since only the 
2^ - 1 subspace is physical. Among the states which span 
the - I Hilbert subspace there are three classes, namely 
all sites carrying the /-particle {ud - 0), all sites carrying 
the Z7o--particle (nd = 1, possible for x > 0) and either b- or 
/- local static occupation (n^ averaged over d'"<->d'^ disorder 
satisfies I > rid > 0). The two latter classes exhibit static spin 
orders or disorders. The appropriate unperturbed ground state 
has to be singled out from these states, and the locally gauge 
invariant Hj mixes within the Qg = 1 subspace all these local 
states, including the d"^«->d^ disorder, in the exact ground and 
excited states. The main formal problem with Eq. Q is that 
the Cu-states are represented in the direct R space while the 
O-states are associated with the reciprocal k space. This re- 
flects the fact that the interaction on the Cu site is large, while 
the other Coulomb^^ or electron-phonon^^ interactions are 
(tacitly) assumed to be small, i.e. that they can be added sub- 
sequently. The result is that the d"'«->d^ disorder is the un- 
avoidable ingredient of the exact many body eigenstates of 
the full slave particle Hamiltonian. This observation^^ will 
be developed here in the fundamental dichotomy between the 
coherent and incoherent correlations. 

Concerning the role of A in Eq. ([T]i it is commonly intro- 
duced as the Lagrange multiplier used to implement the av- 
erage local gauge invariance in the familiar MFSB approxi- 
mation. On the other hand, the time dependent perturbation 
theory in terms of Hj is used here in order to benefit from 
powerful diagrammatic methods. In such theory the constant 
term AN in H plays no role. Such approach is locally gauge 
invariant only asymptotically, i.e. when carried out exactly to 
infinite order In other words, the perturbation sub series leak 
from the Qg-1 subspace and, as will be seen, can even break 
the average local gauge invariance. However it will turn out 
that physically relevant properties are independent of A order 
by order of the Dyson expansions, despite the fact that the av- 
erage local gauge invariance is only approximately satisfied. 

Moreover it will be shown explicitly that these expansions 
converge quickly provided that rtd is sufficiently small, as ex- 
pected on noting that the effects of large Ud are absent when 
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vanishes. Thus, ahhough the theory is formally an ex- 
pansion in terms of tpd, its "small parameter" is a sufficiently 
small nj. The quantitative results, independent of A, are then 
reached in low orders of the Dyson theory. On the other hand, 
for larger n,; it will turn out useful to introduce separate chem- 
ical potentials for the two slave particles. 

It is finally worth to note that, in principle, the time- 
dependent T - Q perturbation theory has advantages over the 
often encountered finite-T Matsubara theory^' This stems 
from the fact that the latter uses the canonical ensemble in the 
full slave particle Hilbert space and thus treats the Qg= 1 and 
the Qg ^ I states on equal footing when they are degenerate 
in energy. It is therefore even more difficult to implement the 
Qg = I local gauge invariance in the Matsubara theory than 
in the time-dependent T - Q approach. One should however 
bear in mind that for low T the Matsubara theory tends to the 
r = perturbation theory and, presumably, most of the low oj 
results derived below for T = may thus be extended to low 
T by the approximate oj <-> iT symmetry. 



2. Unperturbed ground state 

As is well known, the T - Q diagrammatic perturbation 
theory requires that the unperturbed ground state is nonde- 
generate. Keeping in mind that the disordered d'" or d'' states 
are highly degenerate the unperturbed ground state (with total 
spin 5* = 0) has to be chosen either as the n^^ - state or the 
nf* - 1 state with spin order, such as is the Neel state. Leav- 
ing the latter choice for a separate papei'^'we focus here on the 
exact metallic state. We start therefore from the unperturbed, 
tpd - 0, Qj{ - ^, paramagnetic, translationally invariant slave 
particle ground state of Eq. Q equal to 



(3) 



n^?^ = 1 in this state. 
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Consequently, |Gg) is the state with no Zj-particles, i.e. 



JO) 



' - n'j' - 0. Thus rij' = is the outset of our expan- 
sion in terms of small. Since |Gq) is the no-particle state 
it is nondegenerate, irrespective of the Pauli symmetry of the 
feo--particles. 

A side remark is in place here. The above overall state of 
coppers and oxygens is the only state of the Hilbert space with 
required symmetries and has to be used in reaching pertur- 
batively the metallic state with the same symmetries. This 
fact is behind the robustness of our description of the metal- 
lic state. In particular, this state describes straightforwardly 
the CuOi planes doped strongly by electrons or holes, so that 
the magnetic correlations are entirely absent. The magnetic 
fluctuations will be introduced here through the perturbation 
theory on top of the metallic state. This will determine the 
energy scale below which the magnetic fluctuations become 
important. Though magnetic fluctuations are preserving the 
translational invariance in two dimensional systems, for prac- 
tical reasons it is often better to start from the magnetically 



ordered locally gauge invariant states in order to describe the 
behaviors below such energy scale. The resulting coherent 
magnetic fluctuations, such as magnons, restore then a poste- 
riori the translational symmetry. 

Coming back to the translationally invariant state, the /- 
1 Eqs. (2| and (jsjl can be chosen as 



particles in n^"' 



bosons or spinless fermions. The advantage to choose /'s as 
bosons ("slave boson representation", SBR) is that then /jq-'s 
can be taken as fermions indistinguishable from p-fermions. 
This permits to satisfy the anticommutation rules between the 
physical c-fermions on the Cu sites and the p-fermions on the 
O-sites. Such representation is appropriate when starting from 
the Neel state because it can be treated as nondegenerate.^ 
SBR is however disadvantageous for the state (|3]l with bosons 
since the latter is then highly degenerate with respect to mul- 
tiple boson occupations of Cu sites. The difficulty in question 
is obviously eliminated on choosing /'s as spinless fermions 
("slave fermion representation", SFR) and Z^o-'s as bosons, 
since the state of Eq. Q is then nondegenerate. However, 
in that case, we are dealing with three kinds of distinguish- 
able particles, i.e. the anticommutations between the Cu sites 
and O-sites are replaced by commutations. The correspond- 
ing time-dependent slave fermion perturbation theory (SET) 
must be therefore antisymetrized a posteriori (it will be then 
named "antisymmetrized" SET, ASET); this is the route cho- 
sen henceforth. 

An immediate formal benefit of this choice is that the state 
of Eq. ([3]), nondegenerate for spinless fermions, can be sim- 
ply expressed in terms of Eourier transforms /J of the local 

operators /' . Indeed, up to an unimportant phase factor we 
have 



(4) 



where the product over k extends over the whole CuOi Bril- 
louin zone, which corresponds to CuOa unit cells. In other 
words, the Mott state of spinless fermions is equivalent to the 
filled dispersionless band of these particles. 

Similarly, neglecting for simplicity all interactions of p- 
fermions, |Gq) is the usual (nondegenerate) Hartree-Eock 
(HE) state of the p-fermions in the cosine band associated with 
Ep and tpp. This band contains 2«*,"' = 1 + x fermions, asso- 
ciated with the chemical potential /^i+ v. Eor x small this band 
is nearly quarter filled. It is usually folded artificially into the 
Cu02 Brillouin zone in two I, / oxygen bands, anticipating the 
effect of Hi(tpj), which is expected to generate three separate 
bands and make the lowest one nearly half filled. 



3. Full representation in reciprocal space 



Once the unperturbed ground state is expressed^^E] in the 
full momentum representation so should the slave particle 
Hamiltonian - Hqa - AN + Hj. In terms of /J, bi 

k k 

(i - I, I) and their hermitean conjugates we have then 



(Ot 
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(0 (Ot (0 

pk K,(T k,cr 



H, 



. ^ fap!^ + h.c. 

k+q,o- k,o- 



(5) 



V2(|siny| ± Isiny 



assuming the D4 symmetry. Here, tpda^-'\k) describes the fact 
that by annihilating the spinless fermion and by creating the 

bt ^ boson, one annihilates the p^J^ fermion in either of two 

k+q,{r k,CT 

j = I, / bands e^* - Bp± 4\tpp sin (kJ2) sin (ky/2)\. 

The perturbation theory can now be carried out in terms 
of Hj on top of the nondegenerate state of Eqs. (j2]i, (|3]l and 
Q, the normal orderings and the time orderings being well 
defined in the k space together with the Pauli symmetry of 
the relevant b-, /-, and p-particles. Since both the Hamilto- 
nian H^i and the unperturbed ground state are locally gauge 
invariant, translationally invariant on the CuOa lattice, and 
symmetric in spin under time reversal, the SFT will generate 
the exact ground state with the same symmetries, expecting 
that quantum fluctuations prevent the symmetry breakdowns 
in d-2. It is thus left to the SFT to keep = 1, gener- 
ate the LSR nj = n^, obey the anticommutation rules on and 
among the Cu sites, and to satisfy the charge conservation rule 
rid + 2np - I + X. The ASFT is eventually constructed only to 
take care a posteriori of the Cu-O anticommutation rules. As 
will be seen below the procedure will prove independent of A 
in each order. 



4. Elementary propagators 

The elementary bricks which build the time-dependent 
perturbation theory according to Wick's theorem are the 
free-particle propagators. Defining, as usual, t) - 

-iiTb^b'^^it)), we find that the free propagator of the ^-particle 
is dispersionless, 



m 



1 



(jj — Ed — A + irj 



(6) 



Through +iri it describes the intermittent creation of the b- 
particle, while its annihilation is impossible in the no-bosons 
state of Eq. Q. In contrast, the spinless fermions can only be 
annihilated. 



The free propagators of 1 H- jc /^-particles contain both +iri 
and -irj components G'"p^{k, cj) and G^p^(k,oj) according to 
their Fermi distribution in the HF state associated with 

the / - 1,1 bands. The corresponding chemical potential is 
hereafter denoted by f/'^K With x < 1, only the states in the 
/-band are occupied. 

The li-particle propagator is mapped on Dg(t) = 

-(imTZgf^h^gf,(t)bl^(t)}. Df, shown diagramatically 
in Fig.[2j is thus also dispersionless and obtained from Eq. (j6]) 
by replacing Sd + Ahy Sd, 



Df = (co-Sd + iri)-' 



(8) 



thus firstly reproduces nf = = nf and secondly is inde- 
pendent of A, as any physical quantity should be. 



B. Exact pdp and dpd propagators 

1. General expressions 

In the next step we discuss Dy associated with the r-th 
order time-dependent perturbation theory r > 0. According 
to the mapping of the t/-hole on the slave particle-hole pair, 
(lo) is given by the (generalized) Bethe-Salpeter equation 



(r—l) 

Here S_ is the quantity irreducible with respect to cutting 

the /5-lines and (w) is the renormalized four-leg vertex 
given iteratively by the Dyson equation 



rfioj) = rfioj) + rf(u)i:p{aj)rf{co) . (lo) 

in terms of the bare four-leg vertex rl°'((x;) characterized by 

k 

yK^*'*, shown in Fig. [3] 

g2r<"V) = a^Gf<ik, + 4" ' (11) 

j=iJ 

with ffV* given by Eq. (jsj. Eq. (jllj generalizes the bare four- 
leg vertex used previousljEl for tpp = 0. 



Both propagators in Eqs. (|5]l and (|6]l are dispesionless, which 
shows how the approach in the reciprocal space treats the lo- 
cality of the slave particle propagation in the direct space. 



Figure 2: (Color online) Propagator of the spinless fermion 
(green)-boson (blue) pair carrying wave vector k and energy ar, ar- 
rows denote that the /^-particle (blue) can only advance and the /- 
particle (green) recede in time. 
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Actually, t^^^^ in Eq. jioj i can be interpreted as an ap- 
propriately weighted and symmetrized generalization to the 
Emery model of the wide-band propagator on the Ander- 



son lattice. According to Eq. (Ill the intermittent p-particle 



f Jr?' is prepared in two i - 1,1 p-bands, instead of one. 



^pd k is thus the canonical pdp propagator which appears 
naturally (as four-leg vertex) in the perturbation theory for 
the Emery model. For r > 0, this propagation involves both 
/ — I, / p-bands and the degenerate d-state, similarly to Eq. ^ 
for an intermittently created t/-particle. Eq. (fTOb can thus be 



interpreted as the Dyson equation for the f^^r^' ' single particle 



pdp propagator with the Dyson self-energy 



(12) 



irreducible with respect to f^^r^ 



,(0) 



lines. According to Eq. 



for Hj, the lowest order S?' is simply equal to of Fig. |2 

and Eq. (jsjl. In other words D**'' is not only the elementary 
t/-particle propagator but also the essential component of the 
lowest order "local" irreducible self-energy tt^'^ = ^^pd^^^^ ™ 

Eq. (jiojl for f^^rl''. The r - \ procedure thus separates 

out the K-independent free cZ-propagator in the leading 
pc/p-particle self-energy n-*'' on associating in Eq. (Ill the 
A;-dependence of triangular vertices of Fig. [s] with the ap- 
propriate ^-dependent weighting within the free /^-propagator 

Pd k 



issue. Indeed, the SFT removes the d^ state completely, while 
the original large Ud theory treats it as the empty upper Hub- 
bard band at large energy + Ud- Thus, in order to satisfy the 
equality «^'| = n"p_ associated with the full fermion anticom- 
mutation rule on the Cu site and that obtained from D J(t) in 

R 

the f — > 0± limits, one should allow for the additional spectral 
weight at the energy Sd + Ud on the f < side. This adds a 
term 



Jr-l)> 



Li)-{Ed + Ud) + ir] 



(13) 



to S**^ '\ which takes explicitly into account the fact that the 
d^ states are empty. Although permanent (average) occupa- 
tion of the so obtained dispersionless d*^ band is forbidden in 
the f/rf — > oo limit, the full fermion nature of tZ-particles re- 
quires visits of the d** band on the f < side. The point is that 



when the expression ( 13 i is integrated over a> in the Fourier 
transform which determines the overall spectral weight asso- 
ciated with the t < component of S*^''"'' it gives a contri- 
bution fl*''"'* which has to be retained in spite of Ud — » oo. 
In principle, a*'' may be determined from the anticommu- 
tation requirement n^*^ - n^^l of the extended Eq. ^ and is 
expected to be small for r/^^ small. In practice, following the 
spirit of the slave particle theories, we associate the physi- 
cal average occupation n^^ with n^*^ of Eq. (jojl, not worrying 
about cr*''"'*^. This shows however how the SFT handles the 
so called asymptotic freedom. Due to the existence of the in- 
termediate oxygen states such procedure is expected to work 
better for the three-band model than for the large U single 
band Hubbard model. 



Figure 3: (Color online) Four-leg vertex (a)); triangular vertices 



are tpja'''{k) and the red lines are the free propagators Gp (k, w) com- 
bined according to Eq. |TTJ. 



On the other hand (lo) of Eq. (9i describes the cre- 
ation/annihilation of the intermittent c/-particle on the Cu sites 
and its subsequent dpd propagation. The factors S'-f^'^w) in 



(r-l) 



Eq. (9 1 are the same as these involved in F^. Due to this S 

^ J k k 

can be taken as the "effective f ree" pro pagator on the Cu-site 
in analogy with DMFT theories Emm! This allows one to cast 
the Bethe-Salpeter Eq. (|9]l into the Dyson form, in accordance 
with the idea that D^i ^ is the propagator of the single dpd par- 
ticle. 



2. Space-time mapping offermionic properties 



It should be kept in mind that in general D'J' and f ^r'j^ 

k P" k 

do not correspond to fermion propagators, for two reasons. 
The first is that d-d anticommutation on the Cu-site is pro- 
jected on the 2^ = 1 subspace and the second is that d-p 
(anti)commutation rule is not imposed. The latter problem 
will be discussed in Sec. IV Bl while here we turn to the first 



C. Leading correction to single particle p- and d-propagators 

1. SFT construction of the Hartree-Fock propagator 

The last step of the above procedure is to determine the 
chemical potential p'^''^ of the p-fermions in order to satisfy the 
conservation of the total charge 1 +x. In the diagrammatic SFT 
this can be implemented step by step by requiring n*^'^* -H 2n*[* = 
1 -H jc or Irip^ -n^p - x with the numbers of bosons or fermions 
respectively and «^ found below. Alternatively, one can 



(r) 



require n'^ ' + In'p ' - 1 -i- on keeping in mind that the LSR 
Mi - Ud also holds asymptotically in the exact SFT. 

Here we choose the prescription + 2«p ' -\+x because, 
as we shall see now, the r - \ Eqs. Q and ( lOi then become 
equivalent to the HF expression for the free propagators in 
the tpd hybridized Ud - model. The deep reason for this 
far-reaching step is that the Cu site is initially empty, Eqs. ^ 
|4|l, and at the order r = 1 the intermittent particle does not 
probe two-particle effects which involve the Cu site. The sin- 
gle particle propagation is associated with the coherent band 
structure and the band states are subsequently filled on antici- 
pating the local gauge invariance ni, - Ud rather than using the 
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anticommutation rules between the Cu- and O- fermions as in 
the usual HF-theory. 

More precisely, D*.'' involves S?' = D*?' - (oj - Sd + iriy\ 

k k k 

which is associated with unspecified commutation properties 
through the absence of the -irj component. The "single" par- 
ticle problem of anticrossing"*^ between the Ed level and the 
"two" oxygen / = I, / bands is then solved exactly by d'J' and 

rl'* of Eqs. ^ 
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Both t 3r*,'* and in the Dyson form 
exhibit coherent poles belonging to three bands (branches of 



poles) w*/^ denoted respectively by / = L, /, U. The poles u'i 

k k 

in tpd^^ ^6 associated with the residuals (spectral weights) 
z^i\(J{'). They represent thus the spectral weights of the p- 

k k 

particle prepared on O-sites according to Eq. (|TT]i which prop- 
agate in the j-th band. These spectral weights can be ex- 
pressed entirely in terms of the three o}^ . For example, for 
the lowest band L 



U) 



JLd) 



JLp) 



(Sd ■ 



,(LK2 k 



k 



(Ld) 



, (L) (Ox/ (i) (0\ 

(w_ - e ^{oj\ - e i) 

k pk k pk 



(L) 



k k 



k 



(14) 



and similarly for other two bands. According to Eq. (9l z' 
is recognized as the spectral weight of the hole prepared on 
the Cu-site which propagates in the L-band appearing in the 
in D*,'' propagator. Somewhat counter-intuitively zl^' and z'^^* 

k k k 

are proportional. The spectral weight z'i''^^ corresponds to the 

k 

propagator of the hole created symmetrically on two oxygen 
sites. A simple way to determine z^''' is on using the HF sum 



rulezl^'+zl^"' 



1. 



k k 

The chemical potential of the /?-fermions is next de- 
fined as the energy which separates the poles of the pdp- 
propagator f^^rl'' in the upper and lower w-plane. Since 



the same poles appear in D^'', this step fixes their positions 

too, with no reference to the PauU symmetry of the b\f 
pairs. In contrast to which defines the average number 
of p-fermions on the O-sites, /i^'^ can thus be conjugated to 
the average number of particles in the coherent states on O 
and Cu sites. In other words, f/^\Apd,tpd,tpp, x) can be de- 
termined through the approximate charge conservation rule 
n^^' +2n^p^ - l+x, bearing in mind that n*,'^ and n[^'' are defined 
by Eqs. (9 10 1 as functions of the band parameters and 
The whole r - I procedure described above amounts to the 
redistribution of the spectral weights and the Fermi occupa- 
tion factors fi'^ (with accompanying i/T/'s) from two oxygen 

k 

bands i - 1,1 and the empty t/-state into the three coherently 
hybridized bands wV\ j - L,I,U of itinerant states. 



coherently hybridized noninteracting {Ud - 0) pd particles 
with a shift yu"^ -/i**^' in the chemical potential from the upper 
to the lower hybridized hole states. This shift is large when 
single particle anticrossing is important but the variation from 
«^"' = to n^^^ may well be small. 

It is finally worth of noting that the Cu-Cu and Cu-0 anti- 
commutation rules are satisfied automatically in the HF prop- 
agators. In particular, the f — > limit of D^^\t) in the Dyson 

form gives n^^^ - n^J^ - n^J\ i.e. a^^^ = in Eq. Jisl. 



2. Hartree-Fock state and the local gauge invariance 

Although the HF state is obtained here on anticipating the 
local gauge invariance through the LSR Ud - ni, for arbitrary 
n^'' < 1 it will be appear below that the LSR is well satis- 
fied in low order calculations only for n^'' small. Moreover, 
it is immediately evident that local gauge invariance is not 
obeyed exactly for r = 1 (despite 2?' = 1) because the dou- 

R 

ble occupation of the Cu-site is allowed in the pd hybridized 
HF state. However, the singlet HF autocorrelation function 
~ («j V is also small for n'J'' small. It will appear below that 
the LSR and the spin-flip autocorrelation involve the Cu-O an- 
ticommutation in different ways. They will be thus corrected 
in higher orders distinguishing carefully the processes where 
the local gauge invariance replaces the Cu-O anticommutation 
rule from these in which the latter has to be explicitly taken 
into account. 

Since the convergence of the SET is associated in both cases 
with the HF value of n[^'' and since the properties of the HF 
state will prove relevant in higher order calculations too, let 
us now mention briefly here the values of the single particle 
parameters Apd, tpd, and tpp which make n^^\lS.pd,tpd,tpp,x) 
small for a given \ + x. These conditions can be obtained 
straightforwardly on using the spectral weights of D*'' given 
in Eq. ( [T4| . 

The results in question can be conveniently presented on 
giving and dn^J^ /dx for x - 0, having the metallic regime 
for small x in mind. The simplest situatiorP^HlElIIIl corre- 
sponds to tpd » Apd and tpp = 0, where one immediately 
finds n^j^ - 1/2. Note in this respect that /i^'^ coincides for 
X = with the vH en ergy o iyH in the lowest L-band. This 
gives the equal sharing^ESEll of average charge between one 
Cu and two O's. 

Finite tpp can be easily included in this scheme for \tpp\ <sc 
tpd, Apd- The hole doping x - x^h required to reach the vH 
singularity jc,,h was then found analytically'*^ to be propor- 
tional to -32tpp. The large numerical factor multiplying tpp 
in this expression, is due to four tpp bonds per two tpd bonds 
in the CuOi unit cell. Similarly, it can be shown that a rather 
small \tpp\ reduces appreciably the value of n^* for a given x. 

Illustrative results with sizeable \tpp\ can also be obtained 
analytically for 2f^^ = -tppApd, which moreover corresponds 



As announced, the unperturbed ground state of Eq. ^ to the physically relevant band regimd^l. Despite finite f^^ the 



evolves through the prescription «^'^ + 2n]P = 1 + x which intermediate band is then entirely "flat", ajV = Sp, while the 
anticipates the local gauge invariance into the HF state of the two other two bands are simply given by 
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Figure 4: (Color online) The It^^^ = -tppApj example of the flat 
/-band regime in cuprates: dispersions along the high symmetry di- 
rections (Apj = -tpp = V2 tpj, i.e. nj^ = 0.3 at x = 0). The energy 
origin is chosen at and energies are measured in units of V2 tpj. 
The Fermi energy is shown for x = and x = ~ 0.4. 



Figure 5: (Color online) The spectral weights z., on copper site 
plotted along the Fermi surfaces between the antinode (AN) on zone 
main axis and the node (N) on zone diagonal for x = and x = Xy^. 
The band parameters are the same as in Fig. [4] 



k 2 " 2 
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where gj^ are the general^ coherence factors of the 3-band 
model following from Eq. (|5]l. 

This band regime is illustrated in Figs. |4]and|5] It sets ' 
in the L-band with "weak squarelike" Fermi surface'^at the 
energy close below the vH singularities at aiyfj. Noteworthy, 
expanding the square root (when expandable'*^) in the lowest 



band of Eq. 15 the t, t', t" st ructure used in the so called "tight 
binding" band f^iM^M is obtained. Such "single band" 



fits neglect not only two other bands but also replace z^I^^ by 

k 

unity. In contrast the typical behavior of the spectral weights 
of Eq.|l4| along the Fermi surfaces is illustrated in Fig.jsj z^^* 
decreases in the nodal region at the expense of oxygens and 
the projection on oxygens increases with \tpp\. 

The mentioned analytical results can be extended numeri- 
cally to the general values of \tpp\. In particular, n^'' for x = 
is shown in Fig.|6] Similarly, dn^J^ /dx, positive all over the pa- 
rameter space, and the interesting distance w,,// -ju^'^ between 
the van Hove and Fermi energies can also be easily found nu- 
merically. (L)yH - A**'' defines the energy scale Ajp - Sd - 
which will turn out later [see Sec. |V C| in particular] to be quite 
important. 



(16) 



where - w,,// - ([16f^^ + A^^]'^^ - Apd)/2 is independent 

of tpp. Notably, A^p is larger than tpd since at jc = and 
finite tpp falls below cjyff for all choices of band parameters 
[see e.g. Fig.|4j. 

To summarize the results of the above discussion, the range 
of band parameters associated with n^'' < 1/2 at x = is 
determined. As will be argued below, this range defines the 
range of convergence of the present theory, which will be 
shown below to result in > n'^\ It appeared that tpp is quite 
efficient in reducing n^^^ with respect to its tpp - value, due 
to a large numerical factor it caiiies. As already mentioned in 
Sec. II B tpp requires, in contrast to tpd, a finite positive doping 
X - XvH to make the Fermi energy reach the vH singularity. 



IV. PSEUDOPARTICLES GENERATED BY "MIXED 
VALENCE " FLUCTUATIONS 

A. Leading NCA correction to the single slave particle 
propagators and local gauge invariance 

1. Ordinary perturbation theory 

In this Section we turn to the properties of the propagators 



b'"^\1c, oj) and F'"^'{li, lS). As the first they are required to asses 
the accuracy of the above HE result with respect to the local 
gauge requirements = «/, (LSR) and rih + nj - \ (Qg - 1). 
As the second, these propagators will be used to construct the 
next order r = 2 iteration of Eqs. ( |9|10|13| l. 

The disconnected diagrams for b- and f-propagators de- 
scribe the propagation of intermittent /7-particles and /-holes 
as unaffected by the local incoherent p-d "mixed valence" fluc- 
tuations of the vacuum. As usual, these diagrams important 
for understanding the structure of the time dependent pertur- 
bation theory are not shown. The lowest order connected di- 
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Figure 6: ny at x = for general band parameters in particular these 
which correspond to nj,'' < 1/2. 

agrams for AB^P(<jS) and AF'^\(l)) are depicted in Figs. 7 and 
[8] Both diagrams describe respectively how the incoherent 
"mixed valence" fluctuations perturb the propagation of the 
intermittent b-particles and f-holes. The arrows of time, asso- 
ciated with +17] factors in elementary propagators B^^\cij) and 
Ff\(L>) of Eq. ^ are depicted in order to emphasize that tem- 
poral orderings are all important. The bubbles which appear 
in Figs. [7] and [8] are the lowest order irreducible Dyson self 
energies for By(cj) and F^^ \cj). They both contain r*?^(tij) of 
Eq. 
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and involve summation over the occupied states in 
the Z-band as indicated in Figs. |7] and [8] by the /^-propagator 
going (only) backwards in time. 

In Fig. I?] the external frequency enters the bubble ySj,'' from 
right to left, while in Fig. |8]for 0*,'^ it goes from left to right. 
In addition, the spin factor 2 multiplies the /-bubble, unlike 
the ^-bubble, but we remember that there are two fc-bos ons 
for each k. The latter expresses the fundamental property ^3221 
of large Ud situations, namely that each local empty (d'°) site 
can be filled with two (d^) spins. Here large Ud is replaced 
by local gauge invariance, i.e. by the the f ^ b symmetry 
exemplified by Figs. [7] and |8] and the concomitant distribution 
of factors of 2. 




Figure 7: (Color online) Lowest order renormalization for B. The ar- 
rows of time are shown. The energy and momentum enter the bubble 
/j'l" from right to left. 



Figure 8: (Color online) <-» / symmetric lowest order renormaliza- 
tion for F, which defines the bubble <p'j\ 



taken together give An'^p - -An^'\ i.e., AQj^ - 0. However, 
as easily seen, AB^j\a)) and AF^^\(l)) are singular when the 
anticrossin^Sl between the Sd level and bands j - 1,1 affects 
the occupied states in the /-band. Actually, such AZ?||'' and 
AF^^^ coiTespond to the second order Schrodinger perturba- 
tion theory for 3 bands in the "HE propagator" D^'', rather 
then to the exact 3-band theory of Sec. |lllC"T| This leads to a 
large difference between possibly small n^' oif Eq. (|9jl, which 
takes the single particle anticrossing into account exactly, and 
the large number An"^'' of /j-bosons. Such discrepancy cannot 
be removed by taking instead of A«^'^ and the resulting 
convergence of the SET is unsatisfactory. 



2. Perturbation theory re-summed 

When single-particle anticrossing is important it is thus es- 
sential to re-sum the perturbation theory. This can be most 
transparentlyachieved by associating the receding "oxygen" 
line in Figs. Tjandjsjwith hybridized F*.'* rather than with F^"'. 
Concomitantly however, it should be verified that such resum- 
mation is conserving the approximate local gauge invariance. 

Indeed, the /j-bubble is again ^-independent (local), and 
given by 



'_pd_ 
N 



Z 



z'!-'f- 

k ■'k 



CO + (J^^ + A + lirj 



(17) 



where zl^' are the residues of the propagators t ^F*;'* given by 

k k 

Eq. ( [14| i. Eq. ( [17] i is obtained by taking into account the Pauli 
nature of the involved particles, the energy conservation in 
each triangular vertex and by noting that both propagators in 
the /7-bubble are running backwards in time. 

It is instructive to discuss now the LSR = to the 
first order Although n^'' is given by the HE theory and n^'* 
is evaluated in the Appendix it is more insightful to carry their 
comparison as follows. AB*^'* given by Fig. [t] with pdp hy- 



bridized F^'' on the "oxygen" line gives An*,'*, n^'^ can be con- 



veniently determined from the second term 2("*r*,''2*"* of the 
Bethe-Salpeter Eq. Q before its resummation into the Dyson 
form. On noting further that the product of the squared pole 
and a single pole of such 2*;'' occurs also in AB*'\ one finds 



This means indeed that the two diagrams of Figs. |7] and [8] immediately that 
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An 



(1) 



(18) 



Although Eq. (ITSll departs from the LSR rtt, = n,/, the dif- 



ference between An|^'* and n['' is small if n^* is small. As 
discussed in the Appendix, this follows from Eq. ( [T6) l, and the 
observation in Fig. 6 that «[^'' small is equivalent to tpd small 
with respect to A,/^ . The higher order corrections to ^ n 
can then restore the LSR, as the SET converges from Q^^ k 1 
towards the local gauge invariance Qg - I. 

Actually, a similar reasoning applies to the average value of 
Qj^. The spinless fermion self-energy 4>"^^{(jS) has a structure 



If- 

k' k' 



k' -"k' 



■ Sd — A + 2irj 



(19) 



similar to that of -/3*^^\a)) of Eq. |l7jl except for the sign and 
the additional spin factor 2. As there are two bosons per Cu- 



site, the lowest order equation, which is analogous to Eq. (18 1, 
reads again 
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Figure 9: (Color online) The r = I boson and fermion propaga- 
tors Im Bg \ci)) and ImF^Hci)) for x = and for band parameters 



i.e. nj,'' = 0.4 close below nj," = 1/2. (a) 
o(i)> 



^pd/tpd = 3/2, tpp = -tpd, 

denotes the advancing boson pole B^'', (b) the receding pole B^' 
detached from the dense set of receding poles and analogously (c) 
and (d) respectively -Fq"* and -Fq for spinless fermion. In the nu- 
merical calculation the Brillouin zone is sampled by 1024 ^points. 



An 



(1) 



-An 



(1) 



(20) 



The sets of dense poles in Eqs. Q and ([T7| lie respectively 

— (I) (I) 
in the upper (lower) w-half-plane, which makes An^ (^"/ )' 

associated with AB*/' (AF*") of Figs, [t] and [s] finite for 1 +x + 
0, independently of X. 

However, as illustrated in Fig. |9] due to asymmetric roles 
of signs and factors of 2 (boson vs. fermion) in the Dyson 
summation for boson Z?"^- and fermion F*'* -propagators the 
local gauge invariance of Eq. (20i is lost, i.e., 2^'^ 1- This 
is then subject to higher order corrections which asymptoti- 
cally enforce the local gauge invariance Qg-\. Obviously, 
the convergence is fast provided that n*,'' is small. However, 



as shown in the Appendix, on approaching n 



(1) 



» 1 



1/2 from 
n^j\ This is 

a 



below one finds a singular n^ ^ 
closely related to the near coalescence of the advancing and 
receding pole in the boson propagator (e.g. in Fig. [9]l. It is 
then found numerically that n^'* diverges for n^'' « 1/2. As 
was enounced in Sec. llIICTI the latter can be taken to fix the 
radius of absolute convergence of the perturbative SET. The 
extension to larger values of n^'^ requires thus in principle the 
analytical continuation of the exact SET, once summed for n^'* 
small. This however is not attainable in practice. 

Instead, the approximate local gauge invariance 2*!' ^ 1 
(including approximately the LSR) can be implemented in a 
chosen approximation on introducing in Eq. ([TJ the chemical 
potential of slave bosons (additional to X), following a simi- 
lar idea used before in heavy fermion physics.^" In contrast to 
X the boson chemical potential, which enters in the boson site 
energy, separates the poles in the boson propagator which tend 
to coalesce (see Fig. [9]). According to Eq. (A. 6 1, it reduces so 



n*^''. In principle, the boson chemical potential has to evolve 
asymptotically in the exact site energy of the intermittent bo- 
son. The advantage is that the theory is regular order by or- 
der in each order Although useful in practice, such approach 
tends however to conceal the limitations of approximation at 
use. We prefer therefore to continue here with explicit expan- 
sions in terms of small n[^'*. The work with the boson chemical 
potential is thus postponed to the companion paper devoted to 
the quantitative comparison with experiments. 



Beyond these considerations, both propagators B*'' and F^'^ 
turn out here to be dispersionless in the reciprocal space, i.e. 
local in the direct space. It will be argued below that this is a 
general property of slave particle propagators, as suggested a 
long time ago.^' ny and rij, should be therefore understood as 
temporal averages on the given Cu-site. This emphasizes the 
intimate relation between the slave particle locality, causality 
and the local gauge invariance nj + n\-, - 1 . On the other hand 
the LSR relation nj - links the local property nj to n^ 
which involves both the space and time averagings. 



In summary of this Section, the detailed discussion of 
Eqs. ( [T8] l and ( pO| uncovers the mechanism of achieving the 
local gauge invariance in the SET when anticrossing is im- 
portant. The point emphasized here is that quantitatively sig- 
nificant /[-independent results are obtained already in the low 
order SET provided that n^'' is sufliciently small. 
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B. Band narrowing in the physical single particle propagators 
and d"'<-»d'' disorder 



1. Leading NCA self-energy correction 



Once and F\^'((jj) have been determined, they can 

be used to calculate S^'^ ~ (bJ," * F*'^) in Eqs. 1^ and ( 10 1, 
i.e. to advance the iteration one step further to find the prop- 

(2.) (1) 

agators Dl and r_ of the physical particles. Using the re- 

k k 

lation B^j\a>) = Bq\ol> - A) which follows from Eqs. mi and 
( |l7| i, and b ^ f symmetrically F^^\cij) = Fq \'^ ~ we 
find (after integration over u) - X) that S*'^ is independent of 
A. It was pointed out in the preceding section [Sec. IIV A 21 
that B\ {(jS) can be written in terms of one pole in the neg- 
ative w-half-plane and a set of poles in the positive iw-half- 



7(1) 



plane, B*," 



B*,''^ + B?'*^, and f b symmetrically for 



7(1) 



F*''*^ + F^P^ (the superscripts < and > denote arrows 



of time). Importantly, (see Eq. ( |A.2[ )) the relevant contribu- 
tions to the convolution S*'' ~ (By^ * F^'^) come from the 
poles on the opposite sides of the w-axis. 



£(1) ^ iL(b\'Kf\'^) 



1^d(1)> . I7(l)< ^^dX . J7(l)> 



^ -(B\"' * F\ 



FT) 



(21) 



E^'^ is "local" (dispersionless) in the direct space, i.e. it "oc- 
curs" on the Cu site. As shown in the Appendix, we obtain in 
this way, independently of A, 



r(i)> 



(i+«<'V2)4" 



(1) , T 
(jj — + llTj 



(22) 



v(i)< ^ _P± 
A?2 



'■k'.k" 



. . w - w^^i - ul^l +Sd- 4ij] 

k',k" hk' fk" ' 



„(1) 



's, , '^f , c-d, r^i, i„, 0)^'^}, d)^^} in Eqs. (|22|) and (|23|l can all 

o I " bk' fk" 

be expressed in terms of the bare band parameters ; 
X. They are evaluated in the Appendix in the 



1 and doping 
oo limit, 

where the sets of dense poles are treated as cuts in the lo- 
plane, once the temporal decomposition of B*j'', Fj'' and E*'* 
is properly determined. 

The spectral weight of the leading pole E*'^* of Eq. (22 1 
describes the average reduction in the availability of the 
given Cu-site for the free propagation of the additional hole 
due to "mixed valence" fluctuations of the permanent holes. 
1/2 n^'' can thus be understood as the average projec- 
tor which removes the d-states from the coherent propagation. 
The reduction of S^'^* is nominally quadratic in tp^ since such 
are the generic contributions to «*^'* and 1 - nf, as shown by 



Eq. ( A.l l|i. This is accompanied by the shift of the Cu-site 



energy e^y = — jS^^Hoj - e,/ -i- A) where f3\''(Lo - + A) is 



q(i)^ 



negative according to Eq. ( 17 i. 

In contrast to E*'^^, the receding (-irj) continuum of E*''*^ 
describes the dynamic d'°«->d^ disorder. It is shown in the 
Appendix that in the limit — > oo ImE*''^ for n^'' small is a 
step-like function finite in the range 2ju'^''-e^ > at > 2a)M-Sd- 



t*^ in Eq. (23 i for E''^^ is exhibited in order to stress that the 



generic term in the tp^ expansion of E*''^ is nominally quartic 
in tpd which results in the total spectral weight associated with 



(1)9 

the set of dense poles equal to (n , ) [see Eq. ( A.15 i] 



According to the diagrammatic representation of E*''^ that 
contains in particular the convolution of two contributions 
shown in Figs. |7] and [8] this continuum is due to the incoher- 
ent creation and subsequent annihilation of intermediate local 
b^'f -pairs i.e. to the local, dynamic d'"<->d'' disorder, related to 
"mixed valence" fluctuations. It should be noted however that 
the annihilation (creation) of the intermediate boson and the 
creation (annihilation) of the intermediate spinless fermion in 



Eq. ( 23 1, though local, are incompletely correlated within the 
lapse of time between the initial creation of the ^j'Z-pair and 
its final annihilation. E*'^* is thus only the first (NCA) step of 
the perturbation process which eventually ensures that the /- 
particle is always annihilated/created simultaneously with the 
creation/annihilation of the /7-particle. This issue will be taken 
up further at the beginning of Sec.|V A2] 



2. r = 2 pdp propagator 



The "physical" pdp single particle propagator t ^F?' is de- 



termined by E*'* according to Eq. ( 10 1. It will thus undergo 
a modification of the coherent spectral weight and exhibit 
the effects of dynamical d"'<->d^ disorder These effects de- 
scribe how the "mixed valence" fluctuations on Cu-sites, aris- 
ing from the holes permanently in the system, affect the inter- 
mittent hole created or annihilated on O-sites and turn it in a 
pseudoparticle.l^ 



Let us thus consider the pd p-propagator F?' of Eq. (llOb 

in some detail. In the first step we keep E^'*^ but omit E^^. 
Eq. (|l0) for rf combined with E<'> gives then the band nar- 
rowing and the renormalization of the CT gap, while Sp, as 
well as tpp, remain unaffected. 



fW2 

pd 

v(l) 



,(1) 



For small 



= hpd+l3^l\a>^s^}^ +A) 



n*,'V2). Such renormalization 



(24) 



is about half of that predictecPEl by the MFSBT. Concomi- 
tantly, A^^'j is somewhat decreased with respect to Ap^ simi- 
larly to the MFSB theory. 

Next we include E''**^(w) perturbatively and get 
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j,k 



+ z'^^f2,ImS(i)<)2 



where r]s - T]sign(f/^^ 



-PI) 



(25) 



Lj\ 1) and ju^2' remains to be deter- 



mined below, z I Re related to the presumably small [see 
Appendix] (logarithmic) corrections is omitted above for sim- 
plicity. The poles cj i are thus given by the coherent band 



to uF'l are obtained from Eq. upgraded to r = 2. 



narrowing of Eq. (|24|i and the residues ^ which correspond 



ImE*''"^ generates the contribution to f^^r^^' beyond the co- 
herent band narrowing. It should be realized in this respect 
that ImX'^'*^ is a receding term, i.e. carries the same sign 
as ris- This contrasts with the f ^-spinless fermion where 
the appearance of the d'^^-^d^ disorder is related to the dou- 
ble change in sign of the imaginary part of the self-energy 
[see e.g. Fig. [9) and so to the annihilation of /-fermions. 
Note that the sign of ImS*'' is not determined from the fact 
that the step-like continuum falls in Eq. ( [25] l below the Fermi 
level but rather, its receding (-irj) nature follows directly from 
Eq. ( |A.ll| l. This illustrates the importance of the time-space 
analysis, carried out here. 

According to Eqs. ( [T6| and ( [23] l there are several possible 
regimes concerning relative positions of the conduction states 
and the continuum. A particularly simple situatio ns corre- 
As discussed in Sec 



spond to n|^'* small 
can be achieved by the choice of the bare band parameters at 
X = and always by the choice of x < 0. The conditions for 



111 CI 



nV* small 



n^'' small at x 



are shown in Fig. ^ The overall spectral 



weight associated with ImE*''* is then small, equal to (n*,'')^ 



according to Eq. (A. 15 1. Even then ImE^'**^ describes appar- 
ently a qualitatively new stochastic effect. 

In order to illustrate briefly the r - 2 results we choose here 
the case with bare band parameters of Fig. [4] This example is 
interesting for two reasons. First, it shows that the average 
local gauge invariance of the r - 2 term is becoming poor 
already for n^'^ rather small (unless the boson chemical poten- 
tial is introduced). Indeed, the corresponding n^* a; 0.29 for 
jc = is according to the procedure described above accom 
0.59 which differs from 1 



panied by «^'' 



n'-p = 0.16 and 

- 0.36. Second, this example is instructive because in the 
leading nontrivial order the continuum ImS*''^ overlaps with 

(2) 

the bottom of the lowest L-band a>J^ already for r = 2, as 
illustrated in Fig.[TO| 

Apparently, the effect of Im Z*^'^ should be interpreted then 
as inelastic Landau-like damping of the coherent pdp propa- 
gation by incoherent pairs (i.e. by local, dynamic d"'<->d^ 
disorder related to "mixed valence" fluctuations). As well 
known, the Landau damping corresponds to energy, rather 
than to momentum relaxation. This corroborates the ear- 
lier distinction between the pseudo- and quasi- particles in 
cuprates.^ The corresponding pseudo-particle width is then 




Figure 10: (Color online) Coherent and incoherent pdp and dpd sin- 
gle particle spectra stemming from the r = 2 renormalization of the 
HP structure shown in Fig. 4 for nj," « 0.3 at x = 0. In both fig- 
ures the energy origin is chosen at bare and energies are mea- 
sured in units of bare V2 tpj. The leading "mixed valence" effects 
are included respectively by gray (for incoherent dpd-propagator) 
and blue shadings (for pseudo coherent pdp and dpd propagators), 
yu'"', slightly shifted below the /i'*' value, and n^J^ a 0.36 with 
a= 0.08. 



conveniently approximated by a remarkably transparent ex- 
pression 



2n 



f2 iin5:(i)<(^(2)) . 

P'' L,k L,k 



(2) 

Lo The resulting broadening of 

L,k ^ ^ 



(26) 



Im S^''* is taken here at a> 

the Dirac functions into the Lorentzians normalized to unity, 
is shown by the broad line in Fig.[TO]for the pdp propagation 
in the r=2 three-band model. On the other hand according to 
Eq. (jl4jl the weight of the affected d-states is measured by z^^l- 

Turnmg now to the sum rules we note that the Lorentzian 
(2) J L 

^ 13^ (^26 1 does not affect the renormalized HE contribu- 



of Eq. 

tion rip^^'' to np \ For a given chemical potential jl^^ we thus 
have 



„(2) ^ „(2HF) (2mc) 
P P P 



(27) 



where rip'"'^^ is the correction to the Lorentzian approximation. 
This term which includes the variation of ImE*^'*^ aside the 
band can be neglected on the simplest level. 



3. r = 2 dpd-propagator, sum rules and the chemical potential 

In contrast to the pdp propagator, the local, dynami- 
cal d'"<->d'' disorder is entirely revealed in the r=2 dpd- 
propagator given by Bethe-Salpeter Eq. (j9]) which describes 
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the propagation of the b'^ f particle-hole pair. Beside the 
coherent propagation, its first term exhibits ImS^'^*^ of 
Eq. ( |A.ll| i which is finite not only above but also below the 
bottom of the L^^*-band, all over the Brillouin zone in the fre- 
quency range Ium - < co < 2yU^'* - ej, where the coherent 
hybridization is absent. 

Concerning first the coherent dpd propagation Eq. (|9]l 
shows that some of the coherent spectral weight on copper 



is projected out, i.e. depleted to z 



(2d) 



(1 



Z,-. . 



Here z'^^f^^^ is the r = 2 renormaUzed HP spectral weight 



defined by the self-energy S^'^^ of Eq. i22« divided by (1 + 



.(1), 



' l2)n"P and given by the upgraded Eq. 
depletion of the coherent spectral weight z^^"' the first, coher- 



Due to such 



ent term in nf* of Eq. (|28|l is somewhat smaller than n^^ 



ilHF) 



E.g. for n^'^ 



calculated from the renormalized HE bands for a given chem 
ical potential ju*^' 



1 



■ n 



can be used, on linearizing in terms of n 



{2HF) 



(2coh) 



. (2HF) 



(I) 

(rij Y 12. In this sense it can be said that the 
r - 2 HE band plays an auxiliary role. 

In contrast to that, the propagation in the incoherent fre- 
quency range is dominated by the first term Im 2*'^*^ in Eq. (jojl, 
explicitly related to "mixed valence" fluctuations. The local- 
ity of this term is the remnant of the Mott localization within 
the CuOi unit cell, which occurs for other choices of x and/or 
band parameters. In the simplest approximation n9'"'^^ is ob- 



n^'\ to write 



tained by integrating ImE''** of Eqs. (23 i and ( A.ll|i o ver cj. 
This is accomplished in the Appendix [see Eq. (A.15i] with 
the particularly simple result, namely n'^ 



{line) 



which reduces to n 



{line) 



■"''^''^ for n^i* small. 



W 'Y „ 

The coherent S*^'^^ and incoherent X^'^^ components con- 
tribute then essentially additively to average occupations of 



the Cu sites, i.e. n 



(2) 



{IcoH) , {line) 



We have therefore 



(1 + 



)nyn 



+ n 



{line) 



(28) 



Putting these two contributions to together on keeping in 
mind that rij^^^'' depends, as well as np^^\ on the chem- 
ical potential ^^^^ one finds this latter from the sum rule 



,1^ + 2np - \ + X (here again one relies on the local gauge 
invariance rather than on the Cu-O anticommutation rules). 
This is well illustrated for nj^'^ small on writing 



{IHF) 



1 +x 



(2) 



_ . _ (1)2 _ 



(line) 



-4(«^v 



(29) 



f/^'> corresponds thus to the HE chemical potential for the ef- 



fective doping 1 

J2) 



reduced with respect to 1 -H x. Such 
1 + is put into the L*^'-band renormalized according to 



Eq. (|24|. Einally, the value of n?'^^' can be found using the 



renormalized HE results along the lines briefly described in 



Sec. Ill C 1 The conventional Luttinger sum rule for the band 



states (to be distinguished from the LSR = rih) is broken. 
The reason is that the occupied localized disorder states on 
copper carry more spectral weight than transferred on cop- 
per from the occupied itinerant states. The departure from 
the conventional Luttinger sum rule is however small for «^'* 
small. Apparently, the shift f/^^ - which corresponds to 
the shift rfi^^'' — ' is then also small. 

d d 

Importantly, the whole procedure results in n 



(2) 



> n 



(I) 



the local valence fluctuations increase the average copper oc- 
cupation with respect to its HE values n^'^ and n^^^^K 

It can be finally emphasized that the r - 2 results appear 
above as an expansion in terms of n^''. The present theory, 
where d - 2 explicitly, is thus not an expansion in the num- 
ber of dimensions or in large orbital and/or spin degeneracy. 
Even the — > oo limit is unessential here, used only for the 
analytic calculation of the coefficients n>^\ rf^\ e^, ^jlt" ^^^^ 



small Eq. (28 I in Eqs. ([22| and ([23 1. 



4. Iterative slave fermion NCA 



A few general remarks concerning the SET NCA for fi- 



nite but small n^V are appropriate at this point. D^j^' and 



(2) 



r*?' given by Eqs. (122 



23 



25 1 are not sensitive to the omis 
sion of the Cu-O (anti)commutation rules. E*'^ of Eq. (21 



is determined by B^p and F*,'' which in turn are determined 



by yS^j'* and its f b symmetric counterpart tp''^'. The lat- 
ter are determined by F*,'', which is insensitive to the Cu-O 



,(1) 



commutation and, therefore, so are D^^^ and F^'^' provided that 

the relation n^J^ + 2nf^ = I + x, which anticipates the local 
gauge invariance, is implemented. This can be generalized 
to ~ B^p * F^, i.e. the ^ oo NCA is thus insensi- 
tive to the omission of the Cu-O (anti)commutation rules [see 
also Sec. |VB2[ . Moreover, the SET NCA calculation starts 
from the locally gauge invariant state, is / <-> symmet- 
ric for n^'' small, /1-independent and approximately satisfies 
Qj^ ~ I "on average". It can therefore be called a "conserv- 
ing" approximation.'^ 

Assuming that the advancing boson pole (which generates 
the auxiliary HE band) dominates the coherent propagation 
in the r-th order of the SET NCA the latter is expected to 
spread the effects of d"'<->d^ disorder all over the coherent 
spectrum similarly to early slave boson NCA^' (the r = 3 case 
is tractable as will be discussed in the companion paper). E.g. 
for n^'' small Re Y^''^ additional to that of the leading pole can 
be omitted. This amounts to covering the coherent pdp spec- 
trum by blurred tiles of Im E**^' which correspond to distinc- 
tive powers of n^''. The corresponding inelastic Landau-Uke 
damping generalizes straightforwardly Eq. |26) 



-(2) . 



2n 



-=fyUmZ('-"<(w<'-l). 

(<■) L,k L,k 



(30) 



The spectral weight factor of the auxiliary HE band, which 
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appears in Eq. (30i, shows here that the oxygen dominated 
nodal region [see Eq. ([T4]i and Fig. |5) is the least affected by 
the local incoherent "mixed valence" fluctuations. In qualita- 
tive terms this is consistent e.g. with the ARPES result" that 
upon doping the metallic coherence appears most sharply on 
the Fermi arcs around the nodal points. Moreover the damping 
apparently decreases when frequency approaches the Fermi 
level although we cannot determine at present whether or not 
it vanishes there. However, in broad terms the Landau-like 
energy dissipation due to "mixed valence" fluctuations may 
explain the pseudoparticle behavior (also named here the "hid- 
den Fermi liquid" by analogy with the Hubbard model^°). On 
the other hand, the "strange metal" behavior related to the mo- 
mentum (rather than energy) dissipation through various scat- 
tering mechanisms including these behind the "viscous flow 
of charges"— will be discussed elsewhere. 

It is finally worth of noting that, except through sum rules, 
each spin is independent from the other in the SET NCA, as it 
is in large but finite Ud DMET calculations^^* and in the high 
temperature Gutzwiller-like large Ud approaclP^ to d"'<->d^ 
disorder 



5. Comparison to MFSB+fl 

This section is devoted to the comparison of the r - 2 per- 
turbation result and the MFSB-i-fl, which was briefly surveyed 
in Sec. II C 1 and illustrated in Fig. [T| The main point here is 
that, in contrast to the MFSB-i-fl, the appropriately modified 
r - 2 result is smooth through the BR point and that it can be 



reasonably extended beyond n^'' small. 

The r - 2 discussion shows that the regime of small 
corresponds to a positive dn'J' Jdx ^ dn"^^^^ Idx at x small, 
in agreement with Fig. l dn^J^ jdx for x small decreases with 

i 1 ^ 

respect to dny/dx due to correlations, as also found in the 
MFSB+fl result of Fig. [1] In other words, the crossover for x 
small is absent r - 2 and in MFSB+fl for for any combina- 
tion of bare band parameters associated with sufficiently small 



(I) 



n"^\ The coherent r = 2 renormalization gives 

f2^(l + «"V2)ny' somewhat larger then - of the 

MESB theory. On the other hand, the position, the overall 
weight and the structure of the continua differ markedly be- 
tween r = 2 SET and the MFSB+fl of Sec. |ll] [Fig. [T). WeU 
below the Brinkman-Rice (BR) point the MFSB+fl replaces"*^ 
the narrow half-filled "resonant band" above the nan^ow dis- 
persionless continuum, by the partially filled wide band above 
the smeared continuum. However, the MFSB+fl breaks heav- 
ily the symmetry by using large /o - {f + / ' ) (conven- 



tion of Sec. Ill A 1 denotes spinless bosons as /'s all along), 
hybridizes so the /7o--fermions into the coherent HF-states and 
associates the continuum with a single average in the recip- 
rocal space, summing over the occupied b-states. The overall 
weight of the continuum turns so out to be too large. This is 
manifested, before phenomenological correction,**^ in the dras- 
tic breakdown of the LSR = nh- In contrast, the local gauge 
symmetry is properly taken into account by double f ^ b 



are obtained so for n^'^ small, i.e. well below the BR point. 

The other extreme, well above the BR point, corresponds to 
taking «|^'* at x = close to unity i.e. n'p tending to zero in the 
SFT, according to the discussion in the Appendix and Fig. |6] 
Eq. ( |A.l| i then shows that the distance between two poles in 
the B-propagator which moved from the real axis into the up- 
per and lower <y half -plane increases towards infinity but their 
residua stay finite. The contribution of both poles to Z*'' of 
Eq. pT) may be therefore neglected. The coherent Cu-O hop- 
ping is so found to vanish at r = 2. The remaining contin- 
uum shrinks into the receding ("occupied") pseudo-pole close 
below Ed- By Eq. ( |A.15[ ) the overall spectral weight of this 
pseudo-pole is close to unity (for both spins). With hole dop- 
ing yu*^' tends therefore to jump in the ZR manner ''" into the 
oxygen band around Sp for x - Q+. Alternatively, the intro- 
duction of the boson chemical potential separates the coalesc- 
ing boson poles and keeps them on the real axis. The advanc- 
ing boson pole generates then the r - 2 narrow auxiliary HE 
band with spectral weight on copper projected out. The effect 
is similar'*^ to that for resonant band found in MFSB+fl even 
though n^' a: at X = is opposite to the initial value n^"' - 1 . 
Such robustness of our approach can be traced back to the 
fact that the translationally invariant unperturbed ground state 
given by Eq. (|4]) is unique. 

We turn finally to the intermediary case on assuming that 
n^'' at X = approaches 1/2 from below. In the MESB 
the BR point is reached at x = when the boson conden- 
sate (/ + /') 9^ vanishes together with the band-width of 
the conducting band. In the r - 2 SFT the BR point corre- 
sponds roughly, according to the discussion in the Appendix, 
to the coalescence of two poles in the B-propagator and the 
concomitant divergence of the boson number n^''. This in- 
dicates that the higher order perturbation terms are then re- 
quired i.e. that the NCA "mixed valence" fluctuations are en- 
hanced. Such regime can be analyzed on introducing the bo- 
son chemical potential, obtaining so the r = 2 term which is 
finite and the local gauge invariance is approximately obeyed. 
The band-width of the auxiliary HE band, though decreased, 
remains thus finite in the BR point as was anticipated in the 
insert of Fig. [1] 

The rounding up of the BR transition is accompanied by 
the smoothening of other related singular behaviors predicted 
by the MSFT+fl. In particular, dn^J'ldx given by the NCA 
is loosing zeros which appear in the MESB according to 
Such r - 2 results can be employed in the indica- 



Sec. II CI 



tive way close to x^j. For example, one can use n|^'* < 1 /2, 
n^^' < 2/3 estimated at x = for a reasonable guess of bare 
band parameters which correspond to a comparatively small 

rij^^ Idx. 

The smooth regime found here around the BR point 
strongly suggests that Xcs is not related to the (nonmagnetic) 
NCA effects. The relatively strong nearly commensurate mag- 
netic correlations between Cu-sites which set in below x^ are 
expected to increase somewhat the NCA value of i.e. to 

decrease rv^ jdx at Xe,. This indicates that the crossover at 
small Xcs is to be attributed primarily to sizeable AF coher- 
ences. This will be further discussed in the last Sec. I VII 



symmetric r - 2 average of Eq. ( 23 1, and nearly exact results 
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Notably, the r - 2 local continuum lies below the energy 
Sii for arbitrary n'^\ i.e. it contains the contribution of co- 
valent hybridizations within the CUO2 unit cell to the cohe- 
sive energy of cuprates; this helps to explain the profound dif- 
ference between the elastic and plastic properties of cuprates 
and purely ionic crystalPlby placing the former in the regime 
around the BR point, rather than in the ionic limit well above 
it. 

Summarizing, the present r - 2 SFT agrees remotely 
with the picture well known from the MFSB-i-fl [see 
Sec. |11C 1| . Important improvements of the MFSB-nfl are 
however achieved already at the r — 2 level. The renormaliza- 
tion of the coherent single particle propagation differs quali- 
tatively in the two approximations, especially around the "BR 
transition" and serious differences are found in the interpre- 
tation, position, overall weight and frequency structure of the 
incoherent continuum. The SFT cures in part the breakdown 
of the local gauge symmetry in the MFSB-i-fl. The result is 
nearly exact below the BR transition for n^'' small. It allowed 
us to introduce the concept of "mixed valence" fluctuations, 
which appear to be important in general for understanding the 
physics of cuprates. What remains is the inclusion of the mag- 
netic correlations. 



V. MAGNETIC CORRELATIONS IN THE 
ANTISYMMETRIZED SFT (ASFT) 

This section is devoted to the discussion of all vertex cor- 
rections to the SFT NCA with the emphasis on these related to 
magnetic correlations. The problem is again approached from 

(1) 




Figure 1 1 : (Color online) Skeleton contribution (a) to S*^' generated 
by effective two-particle b-f interaction (b). 



p-holes must involve simultaneously a given Cu-site, one by 
the creation of the fe-boson and the other by the simultaneous 
annihilation of the spinless fermion so that the Cu-site car- 
ries either boson or a spinless fermion at any one time. Since 
asynchronous contributions are allowed in the NCA the con 
tribution to 2*^'*^ of Fig 
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the small n , side 



is to be subtracted from Z''* of 
Eq. ( [23] l. This contribution can be easily evaluated noting in 
particular from Eq. (22i that E*^'^* is also nominally propor- 
tional to f^^. However, remembering that the NCA is a con- 
serving approximation for «^'* small (after averaging in time), 
i.e., that the sum rules are satisfied, the contribution of Fig.[TT| 
will be ignored for simplicity in what follows. 



2. Diagonal vertex corrections to E 



A. SFT vertex corrections to single particle propagators 

In this subsection we consider the vertex corrections to the 
self-energy n^''^ which appears in Eq. ( [lo] i for the /j-particle 
propagator. Due to the equality n'-''^ — f-^Z^''"'* of Eq. (12i 
this discussion encompasses also the d-particle propagator o ' 
Eq. (Wb. Since the discussion of the vertex corrections to 0-'''' 
and ^''^ is quite analogous it will be only briefly mentioned 



0<' 
below. 



1. Off-diagonal vertex corrections to E 

The lowest order (skeleton) diagram for S that is not taken 
into account by the NCA is shown in Fig.[TT] As in previous 
low order diagrams, the arrows of time are carefully included 



In this subsection we start by noting that two "squares" Aq 
in Fig. 12 which involve b and / propagators cannot be re- 
duced to the self-energy renormalizations of one of the prop- 
agators appear in the diagram of Fig. [12 



The latter thus extends the present approach beyond the early 



NCA and DMFT results .E'^S'^^^ depends on k, namely it con- 
tains a nonlocal component. In this respect, we emphasize the 
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appearance in Fig. 12 of two particle-particle or particle-hole 
pdpAmts, connecting two f^^A^ 's. In the SFT, the single par- 
ticle lines carry arbitrary spins. The corresponding elementary 
bubbles open singular particle-hole and the particle-particle 
correlation channels. E?* generates the pseudogap effects of 

the magnetic (SDW), charge/bond (CDW/BOW), and Cooper 
(de)pairing correlations. Such terms and their extension to 
higher orders, additional to the r = 1,2,3 terms describing 
the local, dynamic d'*'«-»d'^ disorder, are therefore of qualita- 
tive importance. 

The terms for E^"'^ additional to these shown in Fig. jT2j 



a(0) 



to show the flows of energy throughout the diagram. Fig 
gives straightforwardly a local vertex correction to E^^* related 
to the local effective b-f interaction. It is clear that when 

the anticrossing of p-bands and the <i-state is important the are easily found on noting that the pf^ and B"-^' propagators 

p-propagators in Fig. [TT] have to be interpreted in our itera- rnust occur f b symmetrically in preserving flie lo- 

Such cal gauge invariance. This corresponds to the replacement 
of Aq in Fig. 12 by Aq = AJJ -H Aq (independent of A) for 



2r(i) 



tive scheme as the pdp hybridized propagators f^jF^, 



diagrams take perturbatively into account that two permanent 
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a) 



b) 



Figure 12: (Color online) Skeleton diagrams which contribute to ir- 



reducible E?*: (a) through the particle-particle channel; (b) through 
the particle-hole channel. 




\ / 

V ' V 



Figure 13: (Color online) Skeleton diagrams for two «-> / symmet- 
ric effective p-particle interactions. 



each spin, shown in Fig. 13 The four-leg vertices shown in 
Fig. [13] describe the effective interactions between the pdp 
hybridized particles. Since the ^-dependences are associated 

9 ( 1 ) 

only with the pdp propagating lines f^^r^ , the effective in- 
teraction f^^Ao between the pdp propagators is local. Indeed, 
Aq * are the convolutions of two and two pf* propagators. 



A Ka.b _ J 



led - - W2 



pd 



2/77) 



(31) 



Here, the denote the external frequencies which run coun 



terclockwise around Aq* in Fig. 13 with U] in the upper left 



'4 ^"•'',inde- 



corner, so that oji + 0)2 - oj^ + cda- Remarkably, TjAq 
pendent of A, is the product of four poles. Fig. [13]and Eq. ^ 
show that f^j^Ag* are kinematic interactions, in the sense that 
one particle has to wait for the other to leave a given Cu site 
before it hops to this site. The causal nature of this result is 
associated with the -lirj position of the poles in Eq. pT]|. 

Analogous kinematic terms were first invokecP^^ the 
T — Q theory of transition metals with moderate Ud and 




Figure 14: (Color online) An infinite order expression for S^(<y); the 
first local term corresponds to the NCA. 



extended^^ later to high T^ cuprates. Eq. (31 1 thus represents 
the generalizatiori^^ of this concept to the multiband Emery 
model in the f/^ = 00 limit. 

Actually, the diagrams shown in Fig. [12] completed through 



Fig. 13 are the only skeleton diagrams for based only on 
the vertex f^^Ao. Assuming that the effect of f^^Ao on 



is dominant as is the case for n^'* small (or when enhanced 
by nestinj22lllllll for moderate n^'') it is of some interest to 
construct the corresponding r — > 00 /l-independent expres- 
sion for 2*,°°' on neglecting the off-diagonal b - f interactions 
of Fig. [if] It is also worth to note in this respect that for 

remains finite, in contrast 
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\+ X - Q the vertex f^,^Ao of Fig. 
to to that of Fig.[TT]3. The corresponding infinite order partial 
sum, independent of A, can then be formulated in terms of the 
renormalized value of Aq, A and as follows. 

First, one introduces the infinite order slave-particle Dyson 
propagators in Aq instead of the bare ones, denoting the re- 
sult by Aq. In Fig. [14] this is depicted by the open square. 



Second, as suggested by Fig. 12 and analogously to the usual 



single band perturbation theory, one joins Aq with the fully 
renormalized four-leg vertex A, shaded square, using the pdp- 
bubble formed from two r^Tj^'s, denoted as arrowless in the 

figure. Convoluting the result with f^^F^ one finally obtains 



I.j'icL)), which is depicted symbolically in Fig. 14 It is note 



worthy that the nonlocal nature of then coiTesponds to the 
nonlocal nature of f^^^^^- The backward arrow in the second 
term of Fig. [14] is shown in order to stress that, in general- 
izing the skeleton diagram of Fig. [T2l there must be a pdp- 
line which starts and finishes in Fig. 14 with receding rl"'"^ of 

^ k 

Eq. (Ill, so involving the occupied k's. The second term in 
Fig. 14 thus vanishes for 1 -H x = while the first reduces to 
5^(0) ^ £,(0)^ jjjg should. 

Figure [14] resembles in some respects the exact diagram- 
matic expression for the single particle self-energy in the Hub- 
bard model .^However, this approximation does not take into 
account the Cu-O anticommutation rules for the vertex f^^Ao 
and it takes only partially into account the multisite structure 
of the CUO2 unit cell. The vertex f^,,Ao forbids the simultane- 
ous occupation of the Cu-site by using the spin-caiTying slave 
bosons and spinless fermions while for triplet scattering this 
should follow automatically from the fermionic Cu-O anti- 
commutation. The Cu-O anticommutation issue will be taken 
up in Sec. |VB| while in the following section we continue the 
discussion of the multisite features of the SET in order to grasp 
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better the full structure of the SFT and of the modifications re- 
quired for the inclusion of the Cu-O anticommutation rules. 



3. Diagonal vertex corrections to ji and (f> 

Figure [14 for ^.j^icS) omits the off-diagonal b - f interac- 
tion of Fig. TT] and then it is consistent to omit p - f and 
b - p terms in and 0(w), respectively. However, the 
diagonal b - b "superexchange" interactions stemming from 
Fig. [15] as well as f - f terms, can be included in Fig. [14] 
via the renormalization of single b- and /- particle propa- 
gators analogous to that of Fig. 12 E.g. the self-energy 
is local and these effects, as well as the effects beyond 
Fig. [14] are conveniently discussed using the local compo- 
nent of the transverse spin-flip correlation function t) - 

T I 

Turning thus to we note that the zeroth-order particle- 
hole correlation among the auxiliary spin carriers vanishes 
i.e. xfh = with = of Eq. (jejl. The propagators 



^ gm> + gm< of £q ( [A3} should be therefore used 
instead. The corresponding elementary correlation function, 
proportional to B*'' * B^'', 



2w(l 



is therefore local and linear in n^" ^ to the leading order, 
i.e. negligible for n^'^ small, compared to the corresponding 
correlations x^pl between the /?(i/?-propagators, discussed in 
Sec. 



+ 3iT]co 



(32) 
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Moreover, in sharp contrast with;^^^,'^^, the correlation 



function is dispersionless and incoherent. For oj > the 



continuum in Eq. ([32|) falls in the range e*,V - ojm > co > 



ej^ - ju^'^ and symmetrically for <y < 0, i.e., (w) is real 
around u -0 within the range given by 2e"^^ - yu*''. 

Eq. ( [32| has a transparent meaning. Namely, the auxiliary 
spins on the Cu site cannot rotate freely because hybridized 
with O-sites on the single particle level, x^hl describes then 
the local dynamic spin-flip disorder on the Cu-sites related to 
"mixed valence" fluctuations and is renormalized further, in 
particular by the vertex corrections mentioned above. 

The vertex of Fig. [T4[ which is relevant for local and nonlo- 
— (1) 

cal renormalization of x\,i, has itself local and nonlocal com- 
ponents. In the particle -particle channel the "bare" vertex 
means in particular that two intermittent opposite spin carriers 
on Cu neighbors can hop simultaneously to an intermediate 
O-site, made empty by the removal of two permanent p-holes 
from this site and two spinless fermions from the Cu-sites. 
Apparently the vertex of Fig. [15] generalizes the well known 
superexchange interaction. Notably, it is involving the first 
neighbors if r^F?'* with t„„ = is used, while with r3r*,''*^ 

o pel f, FP pd 1^ 

effective interactions between next-to-next-Cu neighbors are 
generated too. Since the particles on intermediate O-sites are 




Figure 15: (Color online) One of two skeleton diagrams for local and 
nonlocal b — b interactions. The nonlocal component corresponds to 
the superexchange. 



fermions this vertex is not affected by the omission of the Cu- 
O anticommutation rules. 

Since the boson self-energy fi is local it involves only the lo- 
cal component of Xbb- Generally however, the nonlocal renor- 
malization of Xbh should be taken into account before project- 
ing out this component. 

Already this brief discussion allows us to draw two impor- 
tant conclusions. First, the leading spin fluctuations on the 
Cu-sites described by Xbb are local, occur at high frequencies 
and then they are incoherent in time, reflecting the "mixed va- 
lence" fluctuations with spin flip. Second, the vertices which 
describe nonlocal Z?-/? interactions are not affected by the Cu- 
O (anti)commutation rules. Further physically motivated dis- 
cussion of Xbb which includes b - b interactions is postponed 
to the closing Sec.[Vl] while the second conclusion turns us to 
the more formal problem of the Cu-O anticommutation rules. 



B. Antisymmetrized slave fermion theory (ASFT) 

1. Effective repulsion of two pdp-particles 

As we have seen in Sec. |V B 2] the SFT approach generates 
the effective pdp - pdp interaction f^^Ao which, according 
to the above discussion, dominates the coherent correlation 
effects for n^'* small. Its T = structure in time is explicitly 
related to the local gauge invariance. As already mentioned 
several times, its shortcoming is however that it does not take 
properly into account the Cu-O anticommutation rules. On 
the other hand we have seen above that this anticommutation 
is not important for the effective b-b (and / - /) interactions. 

In order to deal with the mentioned problem we remember 
preliminarily that adding intermittently the p- or d-hole in the 
empty band (I + x = 0) leads to the coherent single-particle 
hybridization irrespective of the value of Ud and the anticom- 
mutation rules, as was discussed in connection with Fig. [T4| 
Likewise, one finds then B = B*"' and F = F*^"*, which is anal- 
ogous to the behavior of the phonon propagator in the standard 
polaron problems.'^ 

Next we turn to the two-particle "bipolaron" case. Such ap- 
proach involves the simultaneous creation, on top of the empty 
pdp-hand, of two holes in the singlet or in the triplet state. 
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This is the fundamental problem in which the two require- 
ments of no double occupancy and the d- p anticommutation 
rules come into play simultaneously. Due to spin conservation 



in the triangular vertices of Fig. 13 this interaction turns out 



in the SFT to be the same in the singlet as in the triplet chan- 
nel. As already mentioned, the assistance of slave fermion 
particles is thus required to forbid the double occupation of 
the Cu-site in both spin channels although the Pauli principle 
should do it in the triplet channel. We note further that Aq 
and Aq are both local (dispersionless) in the direct space and 
equal in external frequencies. Though they differ in the sin- 
glet channel by permutation of the outcoming spins they are 
spin independent. In formal terms, the time dependent SFT 
generates the interaction A" + A* which is symmetric with re- 
spect to the permutation of external variables. In contrast to 
that, the perturbation theory for the original Ud Hamiltonian 
(convergent for d=2) has to generate the antisymmetric two- 
hole interaction because^"* two holes with equal spins cannot 
hop simultaneously to the Cu-site due to the Pauli principle, 
irrespective of the values of the band parameters and U4. 

This discrepancy is cured by antisymmetrizing a posteri- 
ori the symmetric SFT scattering vertex AJJ + A^. The anti- 
symmetrization of the spin independent interaction removes 
the triplet scattering, i.e. all scattering diagrams in which 
the two internal /7-lines in vertices of Fig. [T3] carry the same 
spin. On the other hand the singlet scattering stays doubled, 
equal to AJJ + Aq, as it is. This generalizes the result for 
the single-band Hubbard model'^-^or for its continuous free- 
electron limit**^) to the three-band case. The result for two 
holes is readily extended to the multiparticle case because 
the latter can be expressed in terms of the antisymmetric but 
renormalized A (e.g. in Fig. 14 1, again in analogy with the 
Hubbard case^^. The SFT is modified in this way to the ASFT 
which takes care about the Cu-O anticommutation by the an- 
tisymmetrization of the (renormalized) vertices. Remarkably, 
these properties would also follow from the assumption that 
auxiliary /7o--particles anticommute with ;?o--fermions all other 
(anti)commutation rules being kept as they are. We believe 
that this can be generalized to all ASFT diagrams and will 
call the corresponding slave particle representation antisym- 
metrized SFR (ASFR). 



2. Role of antisymmetrization of the pdp - pdp scattering 

Additional insight concerning the role of the antisym- 
metrization is however desirable. Namely, Fig. [T3] can now 
be understood to describe the antisymmetrized (i.e., singlet) 
scattering of two p-particles, which are not completely dis- 
tinguishable because anticommuting between O-sites is taken 
into account from the outset while their single particle prop- 
agation to (and from) Cu-sites is controlled in part by the lo- 
cal gauge invariance rather then by the anticommutation rules. 
The question then arises whether and when precisely the local 
gauge invariance may replace the Cu-O anticommutation in 
the single particle propagation. 

In answering this question we shall first demonstrate for 
r — 2 that essentially the same result is obtained for l}}.^ in 



a) 




b) 



Figure 16: (Color online) Irreducible b <-> / symmetric renormal- 
izations for E'" quadratic in involving: (a) the local renormal- 
ization of the ii-propagator and (b) the local renormalization of the 
/-propagator. 



the ASFT by treating the triplet and singlet vertex corrections 
in the way described above, and for 2^'* in the SFT, where the 
vertex antisymmetrization is not implemented. To this end it is 
important to note that the SFT/ASFT self-energy corrections 
are obtained by closing one of the two p-lines of Figs. [T3^,b, 
as shown in Fig. [2] 



16 



Clearly it corre 



The resulting AS*'' is shown in Fig. 
sponds to the leading contributions to iF' of Eq. ( [2l) given 
by 



AS(i°)> + AI.W< = — (AB<," *F^ + Bl* aM'*) , (33) 



2;r 



where AB 



(1) 



(Bf )Y;\u) and AFf> 



t>^^\cL)) are 

shown in Figs. |7] and [Bland the spin is conserved in each tri- 
angular vertex. Figure[l6| shows that the effective interaction 



Aq* between the hybridized particles occurs even in the lead- 
ing contribution AE*'* to S*''. In the SFT the spin factors of 
the first and the second term of Fig. 



and 
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associated with /3^P(ol>) 
t^^'((jj), are equal to 1 and 2 respectively. The time and 
the Pauli structure of the SFT multiplies the first term by -1. 
The overall result is the subtraction -1+2 = 1 of the two 
relevant terms in Eq. (33 1. In the ASFT the contribution of 



Fig. 16 1 is removed because it involves two bosons of equal 
spin, while the second term [Fig.fTSI^] carries the spin weight 
1 because the contribution with two bosons of equal spin is 
removed from the sum over the spins. This amounts to the 

The overall re- 
AS(i' = ASli', provided that the prefac- 



*'V2 renormalization of the ff^ line. 



suit is again 1, i.e. 
tor of Ay in Fig 



13?, i.e. the singlet scattering in the ASFT 



r ^(Aq + Aq) is itself taken equal to the singlet scattering in the 
SFT. As could be expected, the d-p (anti)commutation rules 
are then entirely UTelevant in AS''' and AS''^ 



M ■ 

Next, the relation Al}^^ - AE<'' is to be extended to the 
Dyson summation which removes the double poles (pf^)^ and 
(B*"')^ from the above discussion. This step is apparently in- 
sensitive to the p-d (anti)commutation rules. In other words, 
by taking l!^ - S*'' the singlet and triplet A^'* are again dis- 
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solved, as in Fig. 16 in the bubble renormalization of B^'^ and 



F*''. At the r = 1, 2 stages of the Dyson perturbation theory it 
is thus irrelevant whether the double occupation of the Cu-site 
with equal spins is forbidden by the local gauge invariance or 
by the Pauli principle. 

In agreement with Sec. |1V B 4[ this reasoning can be gen- 
eralized to the entire NCA subseries because here all interac- 
tion vertices are absorbed in Dyson bubble renormalizations 
of the /-, b- and p-propagators involved. Indeed, by defini- 
tion, the NCA omits all except bubble renormalizations, and 
is therefore the same for th e SFT and the ASFT (but differs 
from slave boson The first term in Fig. 



14 



thus ob- 
viously corresponds to the "exact" SFT/ASFT NCA for the p 



-1) 



(and d) particles through S**^ 

On the other hand the scattering can be explicitly antisym- 
metrized in all correlation processes where it cannot be re- 
duced to the renormalization of the propagators of the slave 
particles. This concerns in particular the second term in 
Fig. 14 Alternatively, following the recipe for the single- 
band Hubbard case,*^^ one can use the unsymmetrized tech- 
nique in this case too, on determining the signs of the so ob- 
tained diagrams by implementing the Cu-O anticommutation 
rule within Wick's theorem. 



proach (in principle supplemented with local vertex correc- 
tions) is comparable to the nonmagnetic^^ DMFT (based on 
single Cu02 unit cell) which assumes large but finite Ud- The 
nonmagnetic spectral broadening appears in both cases and 
the present analysis associates this broadening with dynamic 
d"'<->d^ disorder The ASFT remains a good "conserving" 
approximation**' when it generates nonlocal incommensurate 
magnetic correlations for r > 4. In this respect the ASF T has 
some advantages over the DMFT based on small clusterd^llZZI 
of CuOi units. Such DMFT also corresponds to a selective 
summation of diagrammatic series, as does the ASFT, but 
small clusters are appropriate only for the description of low 
order commensurabilities. The exception is the DMFT ap- 
proach, which allows for general incommensurate magnetic 
correlations'' perturbatively, but it is at present restricted to 
the single band Hubbard model. In the next Section we shall 
therefore examine in some detail the incommensurate, mag- 
netic correlations for the three-band model. 



C. Ladder approximation for coherent magnetic correlations 

1. Elementary ladder summation 



3. Local gauge invariance of the ASFT 

It is desirable in principle to have a locally gauge invari- 
ant theory in which all interaction vertices are properly sym- 
metrized. The interactions between the spin carrying bosons 
or between spinless fermions are antisymmetric already in 
the SFT because mediated by particles on oxygens which are 
fermions. Such is also the case of the interaction between 
the p-particles after the ASFT antisymmetrization. For the 
"off-diagonal" interactions between different particles, such as 
that in Fig. [TT] ^, the proper symmetrization can be achieved, 
if necessary, by the assumption that feo--bosons anticommute 
with pi^-fermions; this ensures the appropriate book-keeping 
of signs within Wick's theorem. However, we have not ex- 
amined this question in detail because it is of limi ted inte rest 
for nj,'^ small enough. Indeed, according to Sec 
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the 

low order ASFT is nearly locally gauge invariant on average 
even with the omission of the off-diagonal vertex corrections. 
Qmr ~ 1 if the f <r^ b symmetry is maintained, step by step, 
in the calculation of the /- and b- propagators. 

Noteworthy, the replacement of the SFT by the ASFT 
which satisfies the Cu-O anticommutation is reflected for 
r > 4 in the on-site spectral weight associated with the empty 
d^ state at "infinite" energy + Ud according to Eq. 



The ASFT theory based on the f/^ ^ oo Hamiltonian of 
Eq. ^ generates the "mixed valence" d'°«->d^ disorder begin- 
ning with the r = 2 single particle propagators and the singlet 
pdp particle-particle and coherent particle-hole correlations 



for r = 4, as shown in Fig. 12 Since the effective interaction 



r^Ao turns out to be repulsive for long times its primary ef- 
fect concerns the coherent spin flip electron-hole correlations 
in addition to the d'"<->d'' disorder Although the d'"<->d^ dis- 
order occurs in the single particle propagation in lower order 
than the spin flip processes, the leading r - 2 disorder con- 
tribution falls far from the Fermi level and the similar obser- 
vation applies to the subsequent r = 3 d'"«->d^ disorder term 
which expands further in terms of Fig 



12 5 will thus be 



first investigated on using the free F*''', Z?**^' and ^pjr^'* prop- 
agators. This amounts to the investigation of the magnetic 
spin-flip susceptibility ^5 in the fully coherent limit. The 
relevant corrections to these results related to the "mixed va- 
lence" fluctuations are briefly examined later. 

The spin-flip response to the transverse staggered magnetic 
field coupled on O sites to the p - d hybridized spin flips is 
associated here with the free bubble ;t'^^'(^, lo). ^ characterizes 
the (intracell) symmetry of the applied magnetic field. For u 
small is dominated by the intraband contribution (denoted 



13| . The from now on as x^f))- Through spectral weights it differs in 

dpd or pdp - pdp bub- 



suppression of the triplet scattering between the pc/p-particles 
replaces (/l-independent) 1}^^ by (/l-independent) ll^^ which 

has a non-local and a local component. The local component is further renormalized according to Fig 



Iff 

general from the corresponding dpd - 
bles [see Eq. ([T4|]. 
(1) 



is involved in n 



n,g anticommutation rule, as conceived 
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inEq. (13i 



~ "dR- 

n^^^ is thus affected beyond the r 

,(4) 



4 level of it- 



A" appears explicitly as the effective forwara 



where 
singlet 



eration though the Dyson structure of D^^^' in Eq. IS mam- 
tained. 

It is finally worth to recall that the NCA SFT/ASFT ap- 



interaction-'^ "'^ between the r - I pdp propagators.^^ 



XSDW 



,y(')+/') 



(1) 



(34) 
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that x^pl in Eq. (|36|l shows peaks for (sufficiently) incommen- 
surate wavevectors so that Ud^ plays (only) the role of the 
forward scattering. This occurs for a sufficiently large im- 




Figure 17: (Color online) ^p^AjJ renormalization of the transverse 
magnetic spin flip susceptibility; spin (vertical arrows) is conserved 
in triangular tpj vertices, while the spin flip occurs in magnetic tri- 
angular vertices (black dots) characterized by the symmetry ^ of the 
staggered magnetic field. 



The actual calculations are simple due to the separability 
of poles in the local interaction f^^A^ of Eq. (jsTJ. It is con- 
venient to carry out the internal integrations first over two in- 
ternal frequencies and then over two internal momenta. For a 
given incoming q, u the real part of these integrals can then be 
written as (x^p Udfj- Assuming that Udf^ is small compared 
to the bandwidth and that external to is small with respect to 
- Bd - /i*'' of Eq. 16 the external frequencies of f^,;Ay in 
Eq. (31 1 fall close to the Fermi energy /i*''. This introduces 
as the cut-off energy in It^^A^ i.e. Ud/j emerges as an "in- 



stantaneous" repulsion with a sizeable cut-off frequency Aj^ 



perfect nesting parameter tpp/tpd- Apparently, Eq. (36i leads 
then to strong incommensurate SDW correlations xs dw for x 
close to the doping x^h required to reach the vH singularity. 
Eq. ( [36) l fixes then the energy scale below which the SDW 
correlations become important. T he lead ing incommensurate 
harmonic turns out to be collineai^3IIIl2^nd the related pseu- 
dogap affects the anti-nodal vH points, leaving the nodal arcs 
essentially in the metallic regime. 

In addition, xf( shows a pronounced maximum at ^ = 0. 
This corresponds to the magnetic effect which generalizes the 
conventional Jahn-Teller rule invoked earUer^^ to explain the 
LTO/LTT transition in some lanthanates. Through the choice 
of the symmetry ^ of the magnetic field in Eq. ( [34| , staggered 
or uniform within the CuOa unit cell, one can compare the rel- 
ative stability of the enhanced spin antiferromagnetism**^ and 
ferromagnetism on two oxygen sites within the unit cell. On 
the oth er han d, the discussion of st ates with permanent orbital 
current^^^'^^' possibly relevani^^'^"' for magnetic properties of 
cuprates requires inclusion of the sizeable off-site Coulomb 
interactions. The detailed analysis of ^ = magnetic correla- 
tions will be further earned out in the companion paper, while 
here we proceed by considering the local renormalizations of 



Ladder summation upgraded 



Ud, 



(35) 



The inequality in Eq. ( 35 i is consistent with the assumption 
that Udfi is small. It follows from the fact that for yU*'^ < oj^h 



Eq. ([T6|) gives A,/^ > 2tpd. 

Such reasoning can be extended to the whole ladder series. 
The result amounts to the renormalization of the vertex f/^^ in 
Eq. S**'' 



Ud^ 



UdnXpp 



(36) 



X'll is to be evaluated with spectral weights z'^' of Eq. (|14|, 



unlike wadx^^p- The same "weak coupling" approxima- 
tion in terms of Udp which replaces f^^A^ by Udp takes also at 
the Fermi level the spectral weights z^i"' involved in;^^^. Ac- 



cording to Eq. ( 12 1 they are equal to zl"' A^^/tj^^j. This means 
that antinodes aominate the pdp - pdp bubble in Eq. (|36]l 
not only due to nesting but also through the enhancement of 
spectral weights z?''' around the van Hove points [see e. g. 



Fig. |5). Further, it is worth noting that Eq. ( [36| omits the 
parqueiSS'^ vertex renormalization and in particular the sin- 
glet Umklapp component of Ud^- This is justified provided 



Eqution ( [36] l can be upgraded on going to the r - \ order 
of the perturbation theory for the slave particle propagators. 
Eq. (A. 5 1 is then used in the upgraded calculation of Ud^i [see 
also the discussion below Fig. [14). Of course, the upgrading 
of Udfi is to be accompanied by the upgrading to r = 2 of the 
/?(i/9-propagators subject in x'^pp '■° '-^^^ effective interaction. 

For present purposes it is sufficient to single out the coher- 



ent contribution to U' 



(1) 



U 



(Vjcoh 
dp 



(37) 



with the depletion (projection) factor familiar from Eqs. p2| ) 
In Eq. (|24jl the depletion factor resulted in the 



and ([24 1 

,2 



(1 + n^i^^ /2)n^'P - f^^'j^ band narrowing. Here, the appear- 



ance of f('>^ in f/("™'' 



pd 

i.. describes the depletion of f^^^-interaction 

Udfi which becomes the interaction of two holes on oxygen 
sites via the intermittently empty copper site. As in the case 
of the band narrowing, the physical meaning of this result is 
quite transparent, namely the reduction of Udp is a direct con- 
sequence of the spectral weight transfer from the coherent to 
incoherent components of the single particle propagators for 
two slave particles. 

It is worth noting that our discussion expresses Udp and 
the ensuing magnetic correlations entirely in terms of the 
(renormalized) single particle parameters as expected for the 
Ud - oo limit. This establishes an intimate relation between 
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the ARPES (related itself directly to the NQR) and magnetic 
measurements which will be used in the companion paper to 
check the present theory against the experimental data. 



3. Comparison to weak coupling and MFSB magnetic results 

The distribution of the roles of the propagators and interac- 
tions in Eqs. (34i and (36i differs from the small Ud theory-''* 
for xsDW and from the ansatz^' which is traditionally used 
in the approach to magnetism^' from the MFSB side. E.g. 



in the weak coupling expression analogous to Eq. (36i the 



repulsion Uj is combined with the dpd — dpd bubble x^^d' 
Similarly, the MFSB ansatTl^il^ uses the replacement Uj — > 
Jpd cos(k^ + cos ky) combined with dpd hubhlex^d^- Note here 
the essential role of the negative sign associated with the cos 
term for k = {n,n~\la, which converts the first neighbor at- 
traction Jpd into the local repulsion. This fundamental step 
was questioned in Ref. HI Equations ([34]) and ( [36| resolve 
thus the long standing contr oversial o f what replaces Jpd of 
the Ud ^ oo Mott-AF limit^EEHH in the well developed 
Ud - oo metallic limit. The answer of Eq. ( 36 1 is that the 



quantity in question is the effective repulsion Uj^ = 2t'^,/A 



combined with the pdp — pdp bubbles x'pl characterized by 



spectral weights z^"' A^^/f^^. Although both proportional to 

fp^ t^rf/j and Jpd are essentially different effective interactions, 
the former between the p-particles and the latter between the 
b-particles as will be further seen in the next subsection. 

Summarizing, once the spectral weights and interactions 
Udp correctly assigned, the ASFT can be restricted to the co- 
herent intraband processes which involve only the "instanta- 
neous" coherent repulsion Udfi in the pt/p-channel. Then the 
tpd, tpp, Udp intraband limit of the 3-band model (with "mixed 
valence" bond fluctuations neglected) maps on the widely em- 
ployed single band f, f', f", Ud weak coupling model.E^EZlThe 
latter can be used to extend straightforwardly the results of 
the present section beyond the ladder approximation. This is 
usually achieved by im proving this approximation in the high 
frequency range,^^^!! which can exceptionally resulP^ in the 
marginal Fermi liquid. Alternatively, anharmonic expansions 
of the Landau type in terms of the magnetic order parameter 
are sometimes used to approach low frequencies"'^^ or else, 
one can start then from the unperturbed static SDW ground 
state with broken translational symmetry, as will be further 
described in the last Section and elaborated in the companion 
paper. 



4. Multicomponent structure of the magnetic susceptibility 

The above discussion emphasizes the role of the single par- 
ticle p-d hybridization and introduces the magnetic suscepti- 
bility related to the resp onse to the staggered field on the O- 
sites. According to Sec. 



VAl 



this is justified for small. 
However, when «^'* ^ 1/2 at jc = the range of applicability 



discussed carefully in function of x. 

This is attempted here by considering the spin flip d-d coi- 
relation function Xddil^ related to the response to the trans- 
verse staggered field on the Cu-sites(it is equal to the longitu- 
dinal one), xlid is containing information about the multisite 
magnetic susceptibility similar to that used semi-empirically 
in the t - J model.l^The advantage of X^jd reveals in 

low orders the magnetic effects on Cu sites. Actually, x\[i is 
connected to a simpler cors:e\?A.or x],],- This follows from the 
fact that the spin-flip operator on the Cu-site i' - c|c| maps 
on b^T-by ff^ in the slave particle theories. Applying Wick's 

theorem to x^jd '^^^ identify two subseries combined with 
the spinless fermion correlations. One of them encompasses 
the pdp - pdp correlations discussed in the previous section. 
The other involves x^ijj, ^ encountered also in the self-energy of 



B-propagators discussed in Sec. VAl 



The leading particle-hole correlator given by Eq, 
well as the generalized superexchange vertex of Fig. [T? 
briefly discussed in Sec. |V A 1| This vertex is nonlocal and 
generates the nonlocal renormalizations ofx),/, itself. The ver- 
tex in question is subject to renormalizations, which introduce 
by the same token RKKY and conventional superexchange 
interactions.^^ To this end the internal spinless fermion lines 
in Fig 



32 as 
were 
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are to be taken as f = f <i'< + F^^> {f ^ b 



symmetrically with Eq. (A. 3 1) rather than as f 



The leading correction to the vertex of Fig. [T4| corresponds 
to the convolution of advancing and receding spinless fermion 
propagators with the corresponding pt/p -propagators. It is lin- 

generates 



ear in «^''. With delocalized pc/p-fermions Fig. 
interactions of RKKY type,'23mediated by pdp- pdp particle- 
hole pairs (with their nesting and interacting properties) which 
couple the spin flips on distant Cu sites. Obviously, this term 
reflect the tuning in of the spin configurations on the Cu and O 
sites additional to that obtained through Ud^i from the simple 
pdp — pdp hybridizations. 

The conventional superexchange interaction J„ d of the per- 
manent (i.e., time averaged) spin^^HolllHllllllll 
neighbor Cu-sites are associated with the transfer of both spins 
to the intermediate empty oxygen site. Both spinless fermions 



15 



in Fig. 15 are then advancing in time. Jpd which appears (to- 
gether with terms which involve more distant neighbors) in- 
volves thus the internal square which is proportional to n 



(1)2 



[see Eq. (A. 11 and below]. In other words it is small for 



(I) 



of the p-d hybridized theory associated with Fig. 14 has to be 



small with respect to the effective interaction f^^Ao between 
the pc/p-propagators (Not surprisingly, just the contrary was 
suggested*''' upon coming from the ti^'' - 1 side). Indeed, the 
two interactions are comparable for n^J^ = 1/2 (i.e. relatively 
large, as noted before^^ for Jpd). However, in sharp contrast 
to fp^Ao, the superexchange interaction involves the F*''* dis- 
order, f ^ b symetrically with the B*''"^ disorder involved in 

In the well developed metallic phase (overdoped and opti- 
mally doped cuprates especially in the high frequency regime) 
the interaction effects on x\\,{q,CL)) come thus out to be small 
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Figure 18: (Color online) Proposed unperturbed ground states in 
function of doping x: (MTL), metallic state, all charges on oxy- 
gens; (N) Neel x = Q state, unity charge on copper, two oxygens 
empty; (ICM's), incommensurate state with AF-magnon frozen on 
unity charge on copper and x on oxygens; ("Neel"), Neel state with 
static d'^^-^d' charge disorder associated with x < 0. Proposed slave 
particle representations for T = diagrammatic perturbation theo- 
ries: (A) Antisymmetrized SFR (ASFR) on top of the MTL state; 
(B) ASFR on top of ICM-states(dominated by {7^.^); (C) ASFR on 
top of ICM states (with U^^ and Jp^ comparable); (D) ASFR or SBR 
on top of the ICM states; (E) slave particle diagrammatic theory in- 
appropriate. 



and incoherent with respect to of the effect of interaction r^Ao 
betwee n the pdp particle-hole propagators discussed in detail 
Notably, x\\,i4^ resulting from the convolu- 



VC 



in Sec. 

tion x^iili^) * Jpd{4y * -A'ti ('^) ill the superexchange chan- 
nel can then be interpreted as describing a dilute gas of RVB 



singletS^^'^ij within the uncorrelated spin disorder of Eq. (32i. 
This picture applies presumably to the optimally and over- 
doped doped cuprates at least for sufficiently high frequencies. 
The summary of this paper is thus devoted to a brief descrip- 
tion of what can be expected on the qualitative level when x 
is decreased from the overdoped to the underdoped x > Xcs 
range in the regime n^J'' ^1/2 for x - 0. 



VI. FROM THE PSEUDO-METAL TO THE MOTT-AF 
INSULATOR 

This section extrapolates in several respects the low order 
nonmagnetic and magnetic ASFT results, found well below 
the BR transition, to its vicinity and is therefore necessarily 
somewhat speculative. Beside the imperfect nesting pa- 
rameter tpp/tpd plays a prominent role in this discussion. The 



results are given in Fig. 18 which summarizes the theoretical 
approach to cuprates proposed here. Note that by w «-> iT 
symmetry Fig.[T8]resembles the well known phase diagram of 
cuprates."* 

We have seen above that elementary interactions, such as 



Udft of Eq. ( 35 1 and Jpd, though nominally comparable in 



value, are quite different in coherence properties when n[^'' at 
x=0 approaches 1 /2 from below, i.e., close to the BR transi- 



correlations are absent and the incoherent "mixed valence" 
fluctuations are the only signature of large Ud. These fluctua- 
tions are localized within the CU-O2 unit cell after excursions 
of the pd/p-particle all over the crystal, as may be expected for 
the nonmagnetic remnants of the Mott phase. This range is 
described here by the Landau-like damping of the single par- 
ticle propagation described within the NCA. Such damping 
gives rise to the pseudoparticle behavioi''". Notably, no con- 
tribution to the resistivity or viscosity^ is expected at low T 
from the energy (rather then momentum) relaxation involved 
in the Landau damping. 

In optimally doped samples |ju - is small, i.e. xsdw 
of Eq. ( [34| becomes large at <d small and its singular behav- 
ior dominates then, through p-d hybridization, the coherent 
magnetic behavior of the system at low frequencies. The vH 
"imperfect nesting" effects tppjtpd, which produce a pseudo- 
gap in the antinodal points, extend thus the doping range of 



magnetic fluctuations in Fig. 18 to relatively large dopings 
X - XyH > Xcs- Landau-like functional expansions can 
be used then to enter the incommensurate magnetic regime 
dominated by Udp- 

When \a)yH - /"I becomes appreciable by varying x the vH 
"nesting" effects in x^pp(4^ became weaker. The resulting 
magnetic dome (B) in Fig. 18 is markedly asymmetric, as 
observed."* Xdd then begins to be driven not only by the hy- 
bridization effect but possibly also through the RKKY inter- 
actions of the Cu-spins involved in Xbh- 

When X reaches Xcs from above ;t'*^|^'(w) and Jpdiq, to) start to 
play a role comparable to that of^p^. Indeed, the incoheren- 
ces of Xhb are related to the dispersion of the L-band. The 
magnetic pseudogap which is a cumulative effect of "coher- 
ent" magnetic correlations, these induced by Udp and by Jpd, 
is likely to progressively localize the pdp-states within CuOa 
unit cells. In particular, the intercell tpp hybridizations are 
then also renormalized down'^i.e. the nodal arcs are flattened 
and gapped away. Consistently, this may also progressively 
activate the Umklapp processes related to Udfi provided that 
tppjtpd is not too large. 

When on decreasing x in Fig. 18 the pseudogap becomes 
comparable to the underlying bandwidth it is appropriate to 
take the (renormalized) pdp propagators as nearly disper- 
sionless. Concomitantly, the width of the continuum in x^^l 
of Eq. (|32]i, modified to include the pseudogap, starts to 
shrink and the leading "mixed valence" spin-flip continuum of 



Eq. ( 32 1 transforms into a pseudo-pole which moves towards 
low frequencies. RKKY interactions are thus gapped away. 
On the other hand, the incoherent superexchange interactions 
between b-particles acquire pseudo-coherence. The resulting 
magnetic coherence, dominated by Jpd, may be considered to 
favorize the long range and long time commensurate ordering 
of the RVB singlets, with weak incommensurability resulting 
in particular from the competition with intercell covalent hy- 
bridizations of doped holes. 

Apparently, the theory is entering so into the magnetic 



tion. In the overdoped range of Fig. 18 the coherent magnetic 



crossover regime close to Xcs, coming in Fig. 18 from high 
frequencies and dopings. On balancing the Mott-AF precur- 
sors which tend to increase «j against the weak decrease of 
NCA n'lj' [see Fig.jrj one may then expect to approach its sat- 
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uration or shallow minimum drid/dx - 0. In other words, the 
whole picture evolves towards the covalently perturbed Mott- 
AF description for x < x^j which is achieved for x = 0. In this 
description the first additional hole goes either to the states 
closer to Sd or to Sp. Both possibilities are named here the 
Mott-AF insulator, since in both cases the permanent hole is 
localized within the Cu02 unit cell (rather than on the Cu site). 
The detailed calculations in the crossover regime are however 
intricate but are expected to yield then a small magnetically 
controlled Xcs- 

Rather than on extending the metallic ASFT beyond its 
range of absolute (nonmagnetic and magnetic) convergence, 
as above, all low energy magnetic crossovers observed in 
cuprates can be alternatively described by changing the na- 
ture of the unperturbed ground state. This is usual in the weak 
coupling limit, except that here one has to be careful to start 
from the locally gauge invariant unperturbed states. In partic- 
ular, for X > one can use as unperturbed states the incom- 
mensurate magnetic states, "frozen magnons" on copper, with 
variable wave vector.ES This is appropriate to describe per- 
turbatively the incommensurate-incommensurate crossover at 
Xcs through the rotation of the leading wave vector from the 
collinear to the diagonal position. Such crossover is expected 
to occur in the underdoped cuprates while the pure metallic 
state is reached in the overdoped range, as shown in Fig. 18 



The case apart is the electron doped {x < 0) magnetic phase 
("Neel"). It contains the (probably quenched) charge disorder 
in the unperturbed ground state which is apparently inconve- 
nient for time dependent perturbation approach with slave par- 
ticles. 

Concerning the choice of slave particle representation, the 
present discussion summarized in Fig. [18] suggests that the 
slave boson representation (SBR) is convenient when the ex- 
act nj is close to unity while antisymmetrized slave fermion 
representation (ASFR) is better when is small. A ques- 
tion open at present is® what representation, ASFR or SBR, 
is preferable at intermediate for x < x^s- The SBR is more 
convenient if the first additional hole goes closer to Sj, while 
the ASFR seems preferable if it goes closer to e,;. This is tan- 
tamount to starting at x = from the AF t - J model, assum- 
ing that the coherent interaction Jp^ (rather than the Umklapp 
Ud/i) plays the most prominent role. For x > x^s the ASFR 
on top of the ICM states might still be better, especially when 
Udft takes over. Noteworthy in this respect, rid in the x > x^s 
incommensurate phase can be increased by the magnetic com- 
mensurability for x < x^j. Such expectation is corroborated 
by the preliminary low order slave particle calculation® for 
X < Xcs which starts perturbatively from the n^"* - 1 Neel 
state and indicates that in the x = "Mott-AF" phase is 
slightly larger than in the small x > x^s metallic phase, 
all parameters of the Emery model except small x being kept 
equal. 

In other words, the crossover at hole doping x^j introduced 
in Sec.|l]on experimental grounds appears here as the result 
of magnetic correlations, separating the insulating magnetic 
phase from the phases with the nodal metallic arc relatively 
weakly affected by doping, i.e. by local "mixed valence" Auc- 



tions. The oxygen dominated nodal arc simply disappears in 
the Mott-AF phase below Xc,. 

The reasoning analogous to that described above for the 
magnetic correlations can be applied to superconducting cor- 
relations too. This however will not be undertaken here, be- 
cause a detailed understanding of magnetic correlations, bet- 
ter than that available at present, is a prerequisite for such an 
analysis. Indeed, in cuprates the high-Tc superconductivity 
is found to occur below the magnetic dome and is confined 
to the vicinity of optimal doping (approximately in range B 
of Fig. [18]) where the coherent magnetic correlations are gov- 
erned, due to the vH (imperfect) nesting, by f/^^. Therefore it 
is not forbidden, as is often advocated in the literature, to think 
that coherent magnetic correlations ("magnons" or "param- 
agnons") are contributing to the glue behind the high-T^ su- 
perconductivity. 

All this deserves further investigations within the approach 
presented here as an expansion in terms of n^J^ small. In 
parallel, the line of reasoning which invokes the effective, 
practically instantaneous local repulsiv e interactio ns between 
covalently hybridized single particles^'"'^'-''!'-'^^^ turned by 
"mixed valence" fluctuations in pseudo-particles'*'^ ^' can be 
checked against the observed properties of cuprates. The de- 
tailed comparison with ARPES, NQR, soft X-ray absorption, 
XPS, x-ray and (non) polarized neutron scattering, Raman and 
optic data including the dc conductivity will be therefore car- 
ried out in the companion paper. It will be argued there on us- 
ing the slave boson chemical potential that the present results 
extended to moderate copper occupations often agree quan- 
titatively with the measurements in the optimally doped and 
overdoped cuprates. The distinction will appear in this way 
between cuprates with tpp comparatively large with respect 
to the observed (renormalized) tpd like YBCCO, BSSCO, in- 
cluding to some extent Bi2102 on one hand, and lanthanates 
on the other hand. However, on decreasing the hole doping all 
cuprates are evolving in the Mott-AF insulators ^'^ ^'^ which 
emerge below Xcs from the metallic phase with "mixed va- 
lence" and magnetic correlations. 
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Appendix: Evaluation of fij,", F,'* and 



(1) 



tuations of Sec. IV B and also by metallic magnetic correla- 



The purpose of this Appendix is to calculate explicitly , 
f\^'> and 2(1) ~ Bf * = b[^^ * F^^ defined in Eq. 
is given by 
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= (A.l) 



CJ-Sd- 



where to' ^ to- A and V) = ySo V') is defined by Eq. ( [l7| . 
Here h-Zt; keeps the memory of the pole in Bf\to) - '(w') 
while the poles in ySg ^(w') are on the opposite side of the to- 
axis. In the usual procedure, one uses the well known equality 



— = P- + indix), 

X + IT] X 



and the oo limit to calculate jSq '(w') in integral form. It 
is th en im portant to retain +iT] in Eq. ( A.l i because Eqs. ( 17 i 
and (A.2i show that lm/3[^\to') vanishes strictly everywhere 



except on a finite segment of the w-axis. ImBj, turns out 
to be negative on the whole Rew'-axis as appropriate for 
bosons. The information about the side of the w' -plane taken 
by the poles is, however, lost by using Eq. (|A.2|i. Analogously, 



Im changes sign twice which is far from the conventional 
(Fermi liquid) behavior of the fermion propagator Additional 
care in the operation with slave particle propagators which 
keeps trace of the pole positions in the complex w' -plane is 
thus required. 

To this end we consider Eq. ( |A.1[ ) as an equation involving 

a discrete set of poles, however dense. According to Eq. ( [T7] i 

(I) — 
the propagator Bl (to ) can be written in the form 



B'il\to') = 



2«+i(w') 



(A.3) 



where « = A^(l H- x)/16 is the number of different poles in f} 



<i) 



of Eq. (A.l I, taking into account that the star of wave vectors 
is degenerate in energy. Then, 



k 



where 1 -i- n'P jl and yf' (independent of id') are the resid- 

" k 

uals of spectral weights associated with the s hifted poles of 
^(dLt'). Indeed, each additive term in Eq. \k.5\ behaves 
manifestly as IJoj' in the \lo'\ — > oo limit, as does the over- 
all expression (A.3i, and therefore the Bq' representation of 
Eq. (A.3 I and Eq. (A. 5 1 are equivalent for our purposes. 

The advantage of Eq. \A.5\ is that the relevant spectral 
weights can be analytically evaluated. In particular, the unique 
advancing pole Bq is well separated from poles in B^' *^ pro- 



(A. 2) is sufliciently small. Its position Sdb is shifted from ac- 



vided that A^^ of Eq. ( 16 is large with respect to tpd i.e. that 



cording to e^^ 



Sd+P^^\t0' 



) i.e. it shifts towards the 

set of n poles. Expanding yS^ '(w') around e^^' the leading pole 
is supplied with the spectral weight 



l+ni>>/2 



(A.6) 



written in the way convenient for the boson on the f < side. 
This determines the average number nj^'' of boson for both 
spins on the Cu-site, which is small for n^'* small. 

On the other hand, although Eq. ( |A.3[ ) is convenient for 
numerical calculations, it is out of reach to evaluate analyti- 
cally the energy shifts Aaf^\k) and the spectral weights y?* 



for each of n poles in Eq. ( |A.3|l for finite A^. Therefore we turn 
to the — > oo limit of Eq. ( A. 5 i. One starts by using Eq. ( A. 2 1 
in Eq. ([l7|for/3*)"(w') 



(A.7) 



It is worthy of note that Eq. ( A.7 1 exhibits the unshifted po- 
sitions and the spectral weights of n poles. These poles are 
grouped in the energy interval [wjvf,//'^'*] where tOM is the en- 
ergy of the lowest occupied pole in the L-band at k - \n , n} 
and //^'^ is the HE chemical potential defined in Sec. 
assert next that in the N ^ oo limit 



IIIC 



We 



i.e., Pn and Qn+\ are polynomials in w' of the order n and n-i- 1, 



respectively. Here s")} and CJ':'' - u/i'' + Ato'i'' are the shifted 



(L) 



XL) 



positions of Sd and w_ 



Bp *(w') is thus a meromorphic func- 



tion which behaves as 1 /to' for \a)'\ — > oo. 

For later convenience the product representation can be 
transformed into the additive form assuming that the renor- 
malized poles remain close to the real axis. This is the case 
for Adfi - Sd - /v*^'* large, i.e. for n^'' small. The result is 



B{'\aj') 
Bll'^icD') 



1+«1'V2 
/ -(I) , • 



+ B'il'<{tv') 



-y- 



k k 



^^oj'- ajf\k) - liT] 



(A.5) 



Imp^^Kto') 



^ (w' - Sd - ReySj,' + (Im/?[; V))2 



= ImB[)'*''(w) 



(A.8) 



using the expression (|A.7|i for \mp^^^(to'). To show this we 



com pare E q. \A.5\ with Eq. \A.%) . According to Eqs. ( |A.5 



and \A.2\ Im B„ can be expressed as the sum of n 6- 



functions, associated with the values of cJ^\k) shifted to 



Kk) 



o^^\k) + Acof\k), 



ImB['^"=(w') 



--7 

N 4^ 



~(L) 



(k)) 



(A.9) 
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The shifts K(J"l^\k) are of the order of 1 jN in the N ^ oo limit 
excepting the two external poles which in principle can detach 
from the rest of the dense set [see Fig.|9). These shifts have 
to b e tak en into account when defining the — > oo limit in 
Eq. (A. 8 This is achieved by performing the ^-integration for 



RcySy' and ImySy ^ in the denominator of Eq. (A. 8 1 while in- 



terpreting ImySg ^ in the numerator as a sum over dense poles 
at the positions infinitesimally shifted to cj^^\k). Such a pro- 
cedure averages over k the ^-functions in the denominator and 
leaves the infinitesimally shifted 5-functions under the sum in 



the numerator of Eq. ( A.8 1 



The analogous procedure can be carried out for the spin- 
less fermion F*''(<y) - Fq ^cj'), as already mentioned in con- 
nection with Eq. coll. In AQ^'^ associated with AbJj'* and 
AFj'* of Figs. iTjandjs] one sums over two bosons according 
to Eq. ( A.ll[ i, and Ag^'* = 0. However, g^' differs from 
unity due to the asymmetry of the signs and roles of the spin 
factor 2 in ' and f'^\ In particular the receding pole in Fg'* 
shifts from zero to -2/3^ '(w' - s^jj), i-C- its distance from 
the continuum increases, in contrast to the boson case [see 
Fig.|9|. Concomitantly, y^^ differs from y^f' of dense states in 



The comparison of Eq. with so interpreted Eq. Eq. jA.lOl l by the appearance of the factor 2 in the denomina 



allows us to identify the spectral weight -/^^ of dense poles 



given by Eq. ( A. 9 1 as 



(1) 



tor, which is unessential. Crudely then, ^ An^' ^ 1 - n 
holds irrespective of the value of n^' 



4" 



„(i) 



Combining the results for nl" and n?' we observe that close 



(b) _ 



to rf^ = 1/2 one has n^'' » 1 - n^'^ * i.e. that bosons 
rather than fermions determine the radius of convergence of 
the SET, as was announced in Sec. |IV A~2l _ 



Re/3<;'(4^'))2 + (Im^<;'(4^'))2 



for each of dense poles, after letting A(JP(k) — > 



(A. 10) 
When- 



Turning now to E*^'* with of Eq. (A. 5 i and its <-> / 



) in the numerator of Eq. ( A.8 i 



symmetric counterpart F^^ known, one finds from S*'' ~ ** 



ever associated with 6iaj' - 

y?^ is obtained by setting a>' = w'i"^ in the 'W-averaged" de- 
nominator of Eq. ( |A.8| l. It follows from Eq. ( |A.8| l that such 
an ImB*,'*^ vanishes outside the interval [wm, In other 



F*" that 



(A. 12) 



words, RejSg' of Eq. (A.lOltakes care of the position and 
the spectral weight of the receding pole detached from the 



using 0o(O) = 2/3o(sd) < under assumption that they are 
small. Moreover, 



dense set, analogously to Eq. ( A. 6 1 for the advancing pole 
By Also, together with Im^Sg \ it gives the "average" spec- 
tral weight y^'^ of all other dense receding poles [see Fig.|9|. 



Once b'q^ = Bq 

r,m> 



'('>>+Bf,"^ known, «!'^ 

' b 



can be evaluated not 



ADAL) (h) if) 
■'ll' ■'P' 'P 'P' 



XL) 



■ oJy' + Ed ■ 
k" 



•4/77 



(A.13) 



only from By on the /> side but also on the /< side as 



s^j^ appears in S*''^ of Eq. ( 22 1 which corresponds to the con- 



volution of only two b f symmetric poles Bq and 
Equation ( |A.13[ ) provides the closed — > 00 expression for 

Ap,P' _ 

Having identified through Eq. (|A.13|), ImE('><(w') 



7(1)< 



lmBf^(oj')doj' 



on inserting the result (|A.10|) for yy in Eq. (|A.ll 



2 
N 

Jb) 



frf, (A.11) 



Not- 
ing that An^'' of Eq. (I81 corresponds to omitting j6[)''(w') 



in Eq. d23]l in terms of the 3-band model parameters. 



can be evaluated according to Eq. ( A. 1 1 1 via Eq. ( A. 2 1, as well 
as its F*'' counterpart. The w-integration of 2 lmTF'{a>) gives 



in the denominator of ■y^\ it is immediately apparent that 

k 



An['' = n"J', provided that n^'' is small 



(1) 



On the other hand, when n^'' = 1/2 is approached from be- 
low the leading pole B'^^^ and the detached pole in By [see 



(2) 

the incoherent contribution to as 



T^- rrn" ' — o t-"-- " "" independent of Sd- This leads for n*,'* small to 

Fig. f)l approach each other due to the logarithmic behavior r « d 



of ReySy According to Eq. (A. 6 the spectral weights of two 



poles diverge at the critical value of n^'' when the two poles in 
question coalesce. Subsequently they split again, leaving the 
real axis in opposite directions. 



by virtue of the approximate local gauge invariance. 



(A. 14) 



(A. 15) 
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